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SECTION I 
DESCRIPTION 



Automatic Polynomial Approximation (SADOl) 

Entire Memory 

Complete, using 1.7 Precision Floating Point Routine A7 

Program requires k minutes 1^ seconds to read in. Time 

required for each problem varies with choice of Mode, 

Option, and Parity, but is given approximately by 

Time Required = [5i + t*t- (Number of Step II iterations )]sec* 

The approximating polynomial of given degree which best 
fits a given function in a given interval in the Chebyshev 
sense is that polynomial of the given degree which has the 
smallest maximum deviation from the given function in the 
given interval. This program may be used to find, approxi- 
mately, that polynomial of given degree which best fits 
the given function in the interval [-1, l]. The user may 
specify the degree of the approximating polynomial desired 
directly or, in Mode I, by specifying the maximum amount 
of error which may be introduced in Step I, as described 
in Section III below. 

The two modes of operation are: 

Mode I : This mode is used if the given function 
is itself a polynomial of degree < ^9 whose coefficients 
are either given on the data tape, are left in the 
machine from a previous problem, or are computed during 
input by an interlude supplied on the data tape. 

Mode II: This mode is used if the given function 
is not itself a polynomial of degree < ^9. Two 
auxiliary subroutines, one to calculate the function 
and the other to calculate its first derivative, must 
either be supplied on the data tape or left in the 
machine from previous problem <, 

In mode I, the preliminary approximating polynomial is 
obtained by polynomial economization; that is, by subtracting 
multiples of appropriate Chebyshev polynomials from the 
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SECTION II 
PREPARATION OF THE 
DATA TAPE 



given polynomial. In Mode II , the preliminary approximating 
polynomial is obtained by interpolation! that is^ by evaluat- 
ing the given function at selected points in the interval 
[-1, l] and then passing a polynomial through these points. 

In either mode, the routine then uses a modified 
Lagrangian interpolation process to find an approximation 
to the given function in [-1, 1] which is of the same degree 
as the preliminary approximating polynomial but which is 
better in the sense that its extrema are more nearly equal. 
This process is repeated until a polynomial is found whose 
extrema differ from each other by less than 1 ^ *, whereupon 
this is taken as the final approximating polynomial. If 
the user desires to use some value other than l^jfor this 
End of Problem Test, he may change the End of Problem Test 
Constant at A7 address ik from "+1-1" = .01 to any other 
desired value by means of a special "Dummy Problem" as 
described in Section V below. 

Numbers on the data tape; must be punched in A7 floating 
point notation, and interludes and auxiliary subroutines 
must be punched in SADOI notation. 

MODE I 
Problem Parameter 1. No. A problem number which 
identifies the problem but does not enter into the computa- 
tion in any way. 

Problem Parameter 2 . m. If 5 (Problem Parameter k 
below) = "+1+152", then the degree of the preliminary and 
final approximating polynomials will be m.: 

This number and 8 together determine the degree 
k of the final approximating polynomial P k (x) in the 
following manner: Since the reduction in Step I from 
A. (x) to the preliminary approximating polynomial Pv( x ) 
■ is actually accomplished by subtracting multiples of 
successive Chebyshev polynomials T. (x), T. -,(x), ... , 
T, , (x) from A. (x), we can predict before the subtraction 
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just what maximum error will result if the proposed 
subtraction is carried out. Therefore, at each stage 
of the reduction in Step I, the routine "looks ahead" 
to see whether or not the proposed subtraction can 
"be carried out without making k less than m or the 
maximum error greater than 8. If this is possible, 
the proposed subtraction is carried out, and the 
whole procedure is then repeated until it is determined 
that the proposed subtraction would violate one or 
the other of the above conditions, whereupon the 
routine proceeds instead to Step II, For further 
details, see Section III, Mode I. 

If the user wants k to be determined by m alone, 
he should set 5 = "+1+152". 

If, on the other hand, he wants k to "be determined 
by 5 alone, he should set m = "+."• 

Note : If Parity (Problem Parameter 3 "below) 
= "+.", then k may be any integer in the range 
1 < k < k9; "but if Parity = "+2.", then k must he an 
integer in the range 2 < k < ^9. 

Problem Parameter j>. Parity. If the coefficients of 
all the odd powers of x are zero, then running time may 
be cut approximately in half by setting Parity = "+2."$ 
otherwise Parity must be set Equal to "+.". 

If Parity is set equal to "+.", no assumptions 
are made regarding the evenness of A.(x). If Parity 

is set equal to "+2.", the given polynomial A. (x) 

1 2 
is assumed to be even, and the transformation y=x 

is used to replace the original polynomial A. (x) of 

even degree and the original interval of approximation 

[-1,1] by a new polynomial Aw^y), where A ± / 2 (y) = 

-x- p 
A. I (x ) = A. (x), and a new interval of approximation 

[0, 1]. Note that the transformation of A.(x) is 
actually accomplished by deleting the coefficients 
ai (i) , a 3 (i) , a 5 (i) , ..., a i-1 (i) of A^x), so that 
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if we set Parity equal to "+2.", we must also liave 
a.^, = aj 1 ^ = ... = a. , ^ = "+." or incorrect 

results will be obtained. ,, 

Problem Parameter ^o <£* If m (Problem Parameter 2 
above) = "+.", then S, the maximum error permitted for the 
preliminary approximating polynomial, determines the degree 
of the preliminary and final approximating polynomials 
as described above. 

Problem Parameter .5. N. The number of subintervals 
into which each of the intervals [-1, 0] and [0, l] is to be 
divided for the purpose of locating the extrema of the error 
• L> £- function. In most cases, an optimum value for N is k, the 
0,^ ^ „ degree of the final approximating polynomial. 



(** , u " The user should be careful to specify an IT 

f /** W ' sufficiently large to insure that all the extrema 

jr jft h will be found, but not so large that the computation 

r -J time becomes unduly long. Too small a choice for 

yi^ U will cause the routine to hang up on an FF 222 Stop; 

$ XJ *~ for a description of appropriate action in this case> 

see Section IV below. 

Problem Parameter 6. D. If D > 0, only the regular 
printout, consisting of the input data and the final 
answers, will be printed; if D < 0, then intermediate results 
are also printed. 

The regular printout consists of the following 
three parts: 

Part I: Problem Parameters 1 through 6. 
Part II: This part of the regular printout, 
depends on which option is specified by Problem Parameter 
7 below. Specifically we have: 
Option A-l. (l) The coefficients of A. (x) read in 

from the data tape. 
Option A-2. (l) The coefficients of A. (x) read in 
from the data tape. 
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(3) The coefficients of A.(x). 
Option B-l. (l) The coefficients of A. (x) left in the 

machine from the last problem. 
Option B-2. (l) The coefficients of A. (x) left in the 

machine from the last problem. 

(2)2* |a ( r V2 r_1 |. 
x ' r=o+l l r ' ' 

(3) The coefficients of A.(x). 

Option C-l» (l) The coefficients of A. (x) computed by 

the interlude. 
Option C-2. (l) The coefficients of A. (x) computed by 

the interlude. 

(2) I 1 .^ |a fr)/^- 1 ! 

v ' r=j+l ' r ' ' 

(3) The coefficients of A.(x). 
Part III: This part of the regular printout also 

depends on which option is specified? 

Option A-l, B-l, or C-l. (l) The coefficients of ^.(x). 

(2) max |K(x)|, min |ET(x) 
x€[-l,l] X xcf-1,1] 1 

(3)Z* |a ( r) /2 r " 1 i. 
x ' r=k+l ' r ' ' 

Option A-2, B-2/ or C-2o (l) The coefficients of P7(x). 



(2) max |E?(x)|, min |£(x) 
xe[-l,l] xe[-l,ll J 

(3)E J Vj _ |a (r) /2 r " 1 |. 

If Problem Parameter 6 is negative, the regular print-. 

out is obtained exactly as described above, but in addition 

between Parts II and III, the routine' prints out the co- 

efficients of each intermediate error polynomial E^(x) or 

E.(x) (omitting coefficients of odd powers of x 
J 
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if Parity =^-2 ."), followed by a list of the extrema 

(x , E. (x )) or (x , E7 (x )) which it has found for 

this polynomial. 

Problem Parameter 7» The user may choose any one of 
six options in Mode I. These are; : 

Option A-l, "Regular" Option, Problem Parameter 7 consists 
merely of the i+1 coefficients of A. (x), where i < 49 and' 



a.^ ' / 0, in order of ascending powers of x. All co- 



(i) 

efficients must be less than "+5+152". 

The maximum total error in this option is given 

by M.T.E. = max \W. (x) | . 
xe[-l,l] * 

Option A-2. Modified "Regular" Option. Problem Parameter 

7 is as follows: 

(1) Polynomial coefficients in order of ascending 
powers of x. 

(2) "+5+155" 

(3) 5 , a new number which replaces the original 5. 

This option enables the user to speed up the solu- 
tion of a "given problem, at some expense in the accuracy 
obtained, as follows : The polynomial economization 
process described in Section III, Mode I, Step I is 
applied to A. (x) to obtain a preliminary approximating 
polynomial A. (x) of degree j, where j is determined by 

m and 5; then the 6 which was originally read in as 

Problem Parameter 4 is replaced by 5 . Then, instead 

of going on to Step II, the routine applies the polynomial 

economization process to A. (x) to obtain a preliminary 

o 

approximating polynomial P, (x) of degree k, where 

K * 
k is now determined by m and o , The routine now pro- 
ceeds to Step II to obtain the final approximating 
polynomial to A. (x) just as though A. (x) were the 
polynomial originally read in from the data tape. 

The maximum total error in this option is given by 
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M.T.E. - Z 1 . -. _ |a (r) /2 r " 1 | + max |£(x)|. 
: r= J +1 r xe[-l,l] J 

Option B-l. "Coefficients Already in Machine" Option, 

Problem Parameter 7 consists of the single number "+1+153"° 

Caution : If the Option B-l problem follows an 

Option A-2 problem, then the polynomial actually 

approximated by the second problem -will not be the 

A. (x) read in by the first problem but rather the 

polynomial A.(x) described under Option A-2. The 

J 
same would be true if the coefficients of A. (x) had 

been computed by an Option C-2 problem or if an Option 

B-2 problem were to intervene between the problem 

which originally furnished the coefficients of A. (x) 

and the Option B-l problem. 

The maximum total error in this option is given by 

M.T.E. = max |e; (x)|. 
x€[-l,l] 1 

Option B-2. Modified "Coefficients already in Machine" 
option. Problem Parameter 7 is as follows: 

(1) "+1+153" 

(2) "+3+153" 

(3) : ° 3 a new number which replaces the original 5. 

This option has: the same effect as Option A-2, in 
that it first causes A.(x) to be replaced by an 
approximation A.(x) of degree j < i, after which the 
routine proceeds as in Option A-2. 

The maximum total error in this option is given 

by M.T.E. = E K^V^ 1 *"" 1 ! + max |£(x) |. 
r_J+1 r xe[-l,l] J 

Option C-l. "Compute Coefficients" Option. Problem 

Parameter ? .consists of: 

(1) "+2+153" 

(2) "00 866K" 

(3) An -Interlude no more than 100 words long 
which enters A7 and computes, for degree i < ^9, i+1 
coefficients and stores them at A7 addresses 1024, 
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1026, „.., 1024 +21; stores Problem Parameter 8 
at 1024 + 2i +2 and ends with a machine language 
"22 649F" order. This interlude may use any Williams 
Memory locations from 866 to 965 inclusive, any A7 
index registers, and symbolic addresses without 
restriction. 

(4) "26 866N" 

After control is returned to the main routine upon 
conclusion of the interlude, the main routine proceeds 
exactly as in Option A-l except that each Coefficient 

a^ is now obtained from A7 address 1024 + 2r instead 

r (ij 

of from the data tape, the last coefficient a. thus 

being obtained from 1024 + 2i and Problem Parameter 8 

from 1024 + 2i + 2. 

The maximum total error in this option is given by 

M.T.E. = max |E. (x) | . 
xe[-l,l] 

Option C-2. Modified "Compute Coefficients" Option. Problem 
Parameter 7 is as follows: 

(1) "+2+153" 

(2) "00866k 

(3) An interlude as described in Option C-l above 
except that the coefficients of A. (x) at 1024 through 
1024 + 2i must be followed by "+5+153" at 1024 + 2i + 2, 
8 at 1024 + 2i + 4, and Problem Parameter 8 at 

1024 + 2i + 6. 

(4) "26 866N" 

This option has the same effect as Option A-2, in 
that it first causes A. (x) to be replaced by an approxima- 
tion A. (x) of degree j < i, after which the routine 

J 
proceeds as in Option A-2. 

The maximum total error in this option is given 

by M.T.E. = Z£ \aJ T ^/2 Tm,1 \ + max |e!(x)|. 
r ~ J+J - r xe[-l,l] J 

Problem Parameter 8. T. The data tape (or, in the 

case of Options C-l and C-2, the list of coefficients at A7 
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addresses from 1024 on) for all problems except Dummy 
Problems must be terminated with either T = "+1+154 " or 
T _ "+2+154". The only difference between these two is 
that in the latter case the routine encounters an OF stop 
before going on to the next problem, 
MODE II. 

Problem Parameter 1. No. Same as Mode I. 

Problem Parameter 2. m. Degree of preliminary and 



v 



A-*" 



final approximating polynomials. 

Caution: Because of difficulties arising from the * $> 
accumulation of round-off error, values of n > 50 are 
not recommended in Mode II. 

Problem Parameter 5° Parity. Set equal to "+•" for 
Mode II. 

Problem Parameter k. 5. Set equal to "+." for Mode II. 

Problem Parameter 5» Same as Mode I. 

Problem Parameter 6. If D > 0, only irlie regular print- 
out, consisting of the input data and the final answers, will 
be printed! if D < 0, then intermediate results are also 
printed. 

The regular printout consists of two parts: 
Part I: Problem Parameters 1 through 6. 
Part II: (l) The coefficients of l£(x) 

(2) max |f(x)-P?(x)|, 
x€[-l,ll ^ 

min |f(x)-P7 (x)|. 
xe[-l,l] 

(3) Zero. 

If Problem Parameter 6 is negative, the regular 

printout is obtained exactly as described above, but 

in addition between Parts I and II, thE routine prints 

out the coefficients of each intermediate approximating 

polynomial P?J(x), together with a list of the extrema 

(x , F(x ) - P (x)) which it has found in [-1,1] for 
r r J£ j 

error function F(x) - P v (x). 



LOCATION 


ORDER 


NOTES PAGE 1 KA 2 






N 




Clear drum 


3 




00 3K 

OOF 00970F 

009K 




Put A7 Floating Accumulator in 
970-975 


9 
10 




229F 509F From Dummy 

Problem 
26l6F 88F i 


Change End of Problem Test 
Constant at A7 address Ik 






8S1^F 8922F 










8J13F OOF 






13 




26762F OOF 
0016K 




Go to (203) to begin new problem 
Put A7 in 16-625 


16 




(A7) 
00832K 




Set Symbolic Address (jl) 


832 


(Jl) 


OOF OOF 
OOQkOK 




Set Symbolic Address (j) 


8^0 


(J) 


OOF OOF 
00852K 




Set Symbolic Address (Nl) 


852 


(HI) 


OOF OOF 
OO866K 




Set Symbolic Address (n) 


866 


(N) 


OOF OOF 






792 


(Yl) 


00792K 
OO75OK 




Put Yl in 792-831 
BEGIN Program E 


750 


(200) 


J0^30F 50(200) 
26 (Yl) 005632F 
00320F 26999F 


From 
SADOI 


. Record Program A on Drum at 
5632-5952 


753 


(200.5) 


J0^-30F 50(200.5) From 
| SADOI 
1 


Record Program B on Drum at 

5953-6273 

. 
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Problem Parameter , 7« The user may choose either of two - 
options in Mode II. These are: 
Option A, "Regular" Option. 

(1) "+5+153" 

(2) "00 65OK" 

(3) An auxiliary subroutine which enters A7, 
leaves F(f ) s where f = (contents of floating accumulator) 
and F is the function being approximated, in the floating 
accumulator, and ends with a machine language "22 621F" 
order • 

(k) "00 700K" 

(5) An auxiliary subroutine which enters A7, 
leaves F*(f), where f = (contents of floating accumu- 
lator), in the floating accumulator i and ends with a 
machine language "22 501F" order. 

(6) "22 65ON" 

These auxiliary subroutines must share Williams 
Memory locations 6^0-7^9 inclusive. They may use A7 
index register b„ and symbolic addresses without 
restriction. 

The maximum total error in this option is given by 

M.T.E. = max |f(x) - p£(x)|. 
x€[-l,l] * 

Option B. "Auxiliary Subroutines Already in Machine" Option. 

Problem Parameter 7 is: 

(1) "+5+153" 

(2 ) "00 65OK" 

(3) "22 63ON" 

The maximum total error in this option is given by 

M.T.E. = max |f(x) - l£(x)|. 
x€[-l,l] * 

Problem Parameter 8. To Same as Mode I. 



SECTION III. 

DETAILED DESCRIPTION MODE I. 

OF METHOD 



In Mode I, the function is given as a polynomial 

A. (x) = Z 1 . a ^x r . where i < k9 and a. ^ £ 0. 

i v r=0 r ' — 1 
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Step I. The Preliminary Approximating Polynomial , 

P°(x). 

First, the routine subtracts appropriate multiples of 

the Chebyshev polynomials T^Cx), T^fe), . .., \ + i( x ^ 

from A. (x) to obtain the preliminary approximating polynomial 

P°(x) of degree k, where k is determined by Problem Parameters 

k 
2 and k (m and o) as follows: 

The routine first reads in the data for the problem 
and checks to make sure that we do not have m > i or 
la.^Vs 1 "* 1 ! > 6 > in eitner of these. cases the routine merely 
skips to the next problem. However, if both m < i 
|a. i V2 i " 1 | < b > then the S iven polynomial 

A. (x) = a n (i) + a, (i) x + ... + a, (i) x 1 

i u 1 1 

is'-f irst reduced to a polynomial of degree i-1 by the formula 

♦ c-i (1) -i (1) &r> xi 



where T (x)=£ t ^^ is the Chebyshev polynomial of degree i, 
i v=0 v 

Note that .since t^ = 2 1 " 1 , A^x) is of degree at most 

i-1, and that since max A i(x)U = 1_ we can use 

xe[-l,l] o 1 " 1 p 1 " 1 



the approximation A ^(x) ~ A ± (x) in the interval [-1,1] with 
the assurance that the absolute value of the maximum total 
error thus introduced at any point in the interval [-1,1] 

a>> 

will be no greater than | ^ |. Note also that, from the 

2 1 

theory of Chebyshev polynomials, if any polynomial with 

leading coefficient a^ other than a ± ([l) -|&i had been 
subtracted from A. (x) to give A ±-1 (x), then the error 
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$ 



introduced in the interval [-1,1] "by the approximation 

a, (±) 

2 
some point in the interval [-1,1]. 



A. (x) ~ A. (x) would have been greater than \~ ± _ ± | at 




Now, : if m > i-1 or a &\ ^ 
|( I !^_| + | 2 \f> 5, then ve 



2" 
cannot subtract: a.^ 1 . ^^ from A ± _i( x ) witnout 

obtaining an approximating polynomial A i _ 2 (x) which is 
either of degree lower than m or has an error greater 
than 8; therefore, in either of these cases, the routine 
now proceeds to Step II. 

r a i (i) *i-i (1_1) 7 

If, on the other hand, m<i-Iani\| ^ | + | r^ 1 J 

< 8, we may subtract a^i ^"i-g fr ° m A i-1^' t0 

obtain a polynomial of degree i-2: 

, % / \ (i-l) i-l(x) 

V 2 W= A i-iW" a i-i IP* 



& ( . (i) ^(x) , (i) (i) Jd: ) 

= \W - a i gTT" " (a i-i a i 2 ±-i ; 

T i-l( x) 

2l " 2 (i) (i-l) „ (i- 2 ) -» 

How, if ■ > i-2 orll^l+1^2— l+l_ ^5~ l J 

/ a _ a 0-) _ili ) -ll 2 - 1/ > 5 then we cannot subtract a ± _ 2 ^ — r-^ — 

i— 1 i o-' - "" o -* 2 



2T 2" 
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from A.' (x) without obtaining an approximating polynomial 
which is either of degree lower than m or has an error 
greater than 8; therefore, in either of these cases, the 
routine now proceeds to Step II. 

r a x (i) vi {w) 

If, on the other hand, m < i-2 and] | - 1 1 |+ | -r-z 1 

(i-2) _ m 

, a i-2 if..* ^ + + (i-2) i-2(x) „ 

+ I . \\ < 5, we may subtract a. p v ' . y ' from 

2 i-3 J *-* 2 ^ 

A. 2 (x), obtaining a polynomial of degree i-3: 



^W-A^W-^"^ 



a) T 



_A 1 tx;-a i ' ^ a i-l a i „i-l 'i-2 



2' 



t .>> ,., ,.,t. ^ 



r a (1) a (I) i-2 . f a (!>_ (D i-1 ■ ) 
" [a i-2 ■ a i 2 i-l ta i-l a i 2 i-l ' 

t l-2 , T i-2(x) 

2 i-2 2 i-3 

This process continues until finally we reach a j such 

(i-v) 
j a._ v 

that m > i-j or L A I "^T r— I > 5, whereupon the routine 

— ° v=0 ' i-v-1 ' * 

takes k = i»-j and sets the preliminary approximating 

polynomial P?(x) = A. . ., v Clearly, the corresponding 
k irJ \*-J • 

error function E' (x) is in this case a polynomial of 
degree i, E°(x) = A ± (x) - A ± (x) = A ± (x) - p£(x), which 

is actually a linear combination of the Chebyshev polynomials 
of degrees k+1 through i. If it should actually happen 

a (i) T 

a. i . / \ 

that E?(x) = — — . ,^ * , then it can be shown that K:(x) 

must have exactly i+2 extrema of equal magnitude and 
alternating sign in the closed interval [-1,1]; namely 
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one at x = -1, one at x = 1, and one at each of the i -. 
points in the open interval (-1,1) for which the derivative 

T. / x 

1 ^ Xj/ vanishes. Furthermore, we know that this polynomial 
dx 

is best in the sense that if ^f . (x) is any other polynomial 
of degree i with the same leading coefficient as E^x), then 

sup |E°(x)|<sup [1R(x)|. 

X€[-l,l] X G [-l,l] 

If, on the other hand, it happens that E ± (x) is some 
linear combination of two or more Chebyshev polynomials, 
it will, in general, not have the above minimal property, 
and in this case it will be possible for us to find another, 
approximation PT(x) to A. (x) such that the error polynomial 

E^(x) = A. (x) - P?;(x) is better than E?(x) in the pense that 
i i k. . -L 

sup |lj(x)| < sup |E?(x)|. 

X€[-l,l] X€[-l,l] 

How the routine obtains this improved approximation is 
described in Step II below,. 

Step II o The Final Approximation Polynomial, p k (x)- 

In Step II, we compute the coefficients of the polynomial 

., E. / \ ^ 

a *^ x ' from the coefficients of E. (x), solve for the roots 
dx i 

of d H*) i n [-1,1], evaluate E. (x) at each such root and 
dx i 

at points x = -1 and x = 1 to get the extrema of E^(x) in 

[-1,1]. If the number L of such extrema is less than k+2, 

the routine encounters an FF222 stop (see FF Stops below). 

If 'C > k+2, the magnitudes of these extrema are compared 

with each other. If these differ -among themselves by 

less than 1 % (this value can be changed by means of a 

special "Dummy Problem" as described in Section V below,) 
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of *^ 

then we take P. (x) = A. (x) - E. (x) as the desired final 

J£ 1 x 

approximating polynomial PT (x) of degree k to A (x), print 
out this result, and go on to the next problem. 

If the extrema differ among themselves by more than 

°( +1 

1 % , the routine attempts to find a polynomial E. (x) 

which is "better than E. (x) in the sense that 

sup |E°f +1 (x)| < sup |E*(x)|. 
xe[-l,l] x€[-l,l] 

To do this, it rejects the smaller (and hence less important) 

of these extrema until only k+2 of alternating sign remain, 

and then finds an approximating polynomial AE (xj of degree 

k to if. (x) such that when AE (x) is subtracted from E. (x), 

we obtain a new error polynomial E. (x) = E. (x) - AE/ (x). 

E.\ (x) also has k+2 extrema in [-1,1] of alternating sign, 

but corresponds to a better approximating polynomial P, (x) 

since its extrema in [-1,1] are more nearly equal than those 

of E.(x). As a result of this process, therefore, we know 

that we .must have 

sup |Ef +1 (x)| < sup \ffW\i 
xe[-l,l] x€[-l,l] 

and that we should actually have 

sup |E +1 (x)| < sup |E (x) | . 
xe[-l,l] x x€[-l,l] 

There are, however, three problems involved in finding 

E. (x): (l) finding the roots of the derivative polynomial 

— ^ X ' in [-1,1], (2) choosing which of the extrema correspond- 
ing to these roots and the points x = -1 and x = 1 to reject 
if there are more than k+2, and (3) using k+2 points to 
determine a polynomial of degree k. 
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The first of these problems is solved "by computing 

—k at the 2N-kL points defined by x = sin I- ^-1= -1, 

. / (N-l) it\ . ( (N-2) A . /N-l Jt\ 

x = sin^- -^Y^gJ' X = Sin (" ^T~2)> ~°> X = Sin l"¥"2-/ 

and x = sin 5- = 1 (or, in case Parity = "+2.", the N+l 

points defined by x = 0, x = sinl- |j, x =1^ |j, .„., 

x = sin [-—— —J, x = sin — = l)| then, starting at the left, 
we test to see whether or not the derivative changes sign 

between each two successive points. (The points x = sinyj ^-J 

w<=re chosen instead of x = | " because- the former are dis- 

T (x) 
tributed in [-1,1] more nearly like the roots of d k % 

dx 

than the latter). If the derivative has the same sign at 

both endpoints of a particular interval, it is assumed that 
, E. / \ 

no root of '— — -, lies between these two endpoints, so 

dx 

we pass on to the next interval. If it does not have the 

same sign at both endpoints, however, then the routine enters 

a root-finding subroutine which locates at least one root, 

-, E . / \ 

say x , of , * in the given interval o The value of the 
r A dx ^ 

corresponding extremum, E. (x ), is now computed and its 

1 r <*/ X 

sign is compared with that of E. (x .. ;, the last pre- 
viously computed extr.emum. If the signs are the same, 
then the extremum which has the smaller magnitude is dis- : 

carded, but if the signs are different, then E. (x J is 

vt <?( r-1 

added to the list of extrema of E. (x) and E. (x ) replaces 

E. (x 1 ) in temporary storage, ready to be compared with 

^v(x n )« This process clearly insures that the list of 
i v r+1 ^ 

extrema of E. (x) so obtained will be of alternating sign 
in [-1,1], and further, if N is sufficiently large, none 
of the extrema will be missed. Note that the extrema which 
must be found in the limit are the k+2 for which the function 
assumes its largest magnitude. Since these alternate in 
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sign they should be fairly well separated. 

The second problem is solved as follows: Let the list 
of extrema of E. (x) found by the above process consist of 

the values m£ = |E ± (x 1 )|, M^ = ^(x^, ..., M^ = |E ± (x^)|. 

(l) IfA = k+2, the routine proceeds to (2). If C = k+3, 
the routine rejects the smaller of M^ and M^i (if M^ = Ma , 
it arbitrarily rejects MVfl ) and proceeds to (2). 

Ifj}> k+k, the M.'s are arranged into £ pairs as 
follows: There are 11 -1 pairs ^ , v 1 J, V> 2 , v 2 J, ..«., 
(jifl , v f)_i) S iven ^y the formulas /i. = max(M , M +1 ) 

ex' . pi bl . ■ . / ^ «^ \ . 

and v. = min (M., M. +1 ), and one pair,, ()i/j\ , v^\ ) given by 

o( , of o( x c< , <K <X N 

/i£ = max (Ma , M^ ) and V/£ = min (M£ , M^ ). 

The routine now rejects that pair of extrema having 
the smallest v. , except that if two or more of the pairs 
have equal v^'s', it then rejects., from among those with 
equal v. 's, that pair which has the smallest p $ or, if 
two or more of these pairs have both v? and^. equal, the 
routine rejects the leftmost pair| i.e., the pair correspond- 
ing to the smallest values of x. Now, ^ is replaced by 
L -2 and the routine goes back to (l). Clearly, the al- 
ternating sign property of the extrema is preserved by 
the above process, since the extrema are rejected in consecu- 
tive pairs. 

(2) The third problem is now solved as follows: although 
it is not possible in general to pass a polynomial of degree 
k through k+2 points it is possible to determine a number h 
and a polynomial of degree k such that the residual of the 
polynomial with respect to each of the points is equal to 
+h or -h and such that the residuals alternate in sign* 

Accordingly, the routine first determines the number h 
by setting the (k+3) st divided difference of the points 

(x x , E. (x x ) + h), (x 2 , E. (x 2 ) - h), ..., (x^, E_. (x k+2 ) 
+ (-1) + ~4i) equal to 0, which gives 



- 18 - 
k+2 E*(x ) *+2 

v JC ' (X ) i=1 1 

h = _ JL± I— , where *« (x ) = lim — / ) : 

k+2 (-1) r 

rSi *' (x r ) 



(x - x 1 ) ... (x- x r _l^ X " x r+i ) ••• ( x ' x k+2^ 

_ , . . th 
Having found h, it now finds the coefficients of the k — 

degree polynomial which passes through the k+1 points 
(x x , E^) + h), ..., (x k +1 , E. (x k +1 ) + (-l) 3 ^) 

by means of the formula 



AE* (x) =([( Zl (x-x 2 ) + z 2 )(x-x 3 ) + z 5 l (x-x^) + ■ ... + z^ (x-x^) + z^, 



where ^ k 

z o, + i= E i (x k + i )+ ^ ,> 



* -. -•» 



( x fc - x k.+l- ) (x k+l " X fc ; 



E^( Xl ) + h E ^(x 2 )-h E f(*E + l )+ (^ 



It will be noted that each z. is exactly the (k+1- j ) — 

J ex' , >. 

divided difference of the points (x k , +1 _j> E ± ( x k +i-j' + 

(-l) k - J h), ..., (x t+v E^ (x k+1 ) : *'(.l)V: 

and also that this formula is actually equivalent to the 

lagrangian interpolation formula : 

"k+1 



v k+1 • . i (x-x ) 



^(x)-£ =1 (E-(x.) + (-l)^) ^ x) J (x ) 

J J 

but has the advantage that it can be evaluated using only k 
multipli 
formula . 
HaV 
AE?(x), we now compute the coefficients of E^ (x) = E ± (x) 



p 
multiplications instead pf the k required by the Lagrangian 



Having thus obtained the coefficients of the polynomial 

<* +1 * 
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AELCx,} and go back to the "beginning of Step II, with °< 
replaced by its new value r ^\+l. 

MODE II 
In this mode the function to be approximated is 
specified by means of auxiliary subroutines for F(x) and 

~, K ' instead of being a polynomial specified by its 
dx 

coefficients . 

Step I . After the problem parameters and the auxiliary 
subroutines have been read in, we do not subtract off 
Chebyshev polynomials from the given function as in Mode I 
"but" instead set k = m immediately and compute F(x) at the 
k+1 points x x = - sm^^-), x 2 = - sin^^, ... > 

x k - sin (kTi i)> vi - sin (kli £)* then obtain *£>> as 

the polynomial which passes exactly through the k+l(=m+l) 
points (x-^ F(x x ), ..., (x fci+1> F ( x k+ i-))* 

Step II. We now proceed as in Mode I, Step II, except 
that our error function #\x) = F(x) - Fj-.-M is no longer 
necessarily a polynomial and therefore is given at a 
particular value of 2£>as the difference of two functions 

F(x) and P 1 (x) which are evaluated separately. Also, 
k o(. 

oC P (x) / \ 

d ^ ^ is given by subtracting ^—^ from d I r , where 

dx clx ax 

the two derivatives are also evaluated separately. Except 

for this difference in the method of computing E (xj and 

. > J , however, the procedure for Mode II, Step II is 

dx * 

exactly the same as for Mode I, Step II. 



LOCATION 



754 



760 



ORDER 



(201) 



(202) 



(203) 



770 



(204) 



(205) 



780 



(206) 



(207) 



(208) 



(209) 



26 (Yl) 005953F 
00520F 26999F 
J0430F 50(201) 
26 (Yl) 0062 7^F 
00320F 26999F 
J0430F 50(202) 
26 (Yl) 006595F 
00320F 26999F 
50430F 50(203) 

26 (Yl) OO5632F 
00520F 26626F 
50^50F 50(204) 
26 (Yl) 005632F 
OOF OOF 
50^50F 50(205) 
26 (Yl) 005953F 
00320F 26626F 
50430F 50(206) 
26 (Yl) 005953F 
00320F 50(207) 

OOF OOF 
50430F 50(208) 
26 (Yl) 00627^F 
00320F 26430F 
50430F 50(209) 
26 (Yl) 006595F 
00320F 26430F 



From SADOI 



From SADOI 



From SADOI,, 
2 (28.6), (28.8), 
(64. 6) 



From (71.57), 
(72.6) 



By (71.56), 

(26) 

From 2(23) 



NOTES PAGE 2 ■ KA 2 



From 1(71.55). 
(126), 3 (26 .2) 



By (71.55), 
(30.10) 
From (4-2.62) 



Record Program C on Drum at 
627^-6594 



Record Program D on Drum" at 
6595-6915 



Playback Program A 



Go to (7.004) to begin first 

problem 

Playback Program A 



(Overwritten after Read -in of 
Program D) Go to (27. l) 
Playback Program B 



Go to (28.5) for first itera- 
tion of Step II 
Playback Program B 



From (71-02) 



(Overwritten after Read -in of 
Program D) Go to (30.05 ) 
Playback Program C 



Go to 0*2.65) to continue Step II 
Playback Program D 



Go to (71.5) to continue Step 



II 



Section IV 
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STOPS 



2V Stops 

I. Control Panel: 20 272 Program: (7-004) 20 (7.004) 

50 272 50 (7.004) 

273 
This stop signifies that the routine has completed the 
current problem. A Black Switch Start will now cause the 
routine to go on to the next problem o 

II. Control Panel: 24 021 Program: l(210) 24 33? 

26 272 26 626F 

30F 

This stop was included so that the operator could run 

an Option B-l or B-2 problem without re-running all of the 

Option A-lorC-1 problem which originally specified the 

(4) 
coefficients a v '.. Therefore, the stop occurs immediately 
r /. \ 

after all the a ^ L ' have been read in (Option A-l) or 
r 

computed (Option C-l). Jl Black Switch Start will now cause 
the routine to continue the Option A-l or C-l problem, 
whereas repositioning the data tape to the beginning of an 
Option B-l or B-2 problem and White Switching will cause 
the routine to run the Option B-l or B-2 problem using the 
coefficients A.(x). 

III . Control Panel : 24 1J J Program : (71 - I 6 ) 24 (71 • 78 ) 

22 1JJ 22 (71-78) 

1JJ 

This stop was included so that the routine could be 

made to execute extra Step II iterations even though the 

End of Problem Test has been satisfied. 

If the extrema do not satisfy the End of Problem Test, 

no stop occurs, and the routine automatically goes on to 

0^ +1 
compute a new error function E (x). If the End of Problem 

Test is met, however, a "Black Switch Start at this point 

will cause the routine to finish the current problem and 

go on to the next- A White Switch Start, however, causes 

the routine to execute an extrai Step II iteration followed 

by another 24 1JJ stop. 
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The routine may be made to execute as many extra itera- 
tions as desired by White Switching once for each extra 
iteration . 

OF Stops 
I. Control Panel: OF 000 Program: (28) OF F 

26 272 26 (7.00^) 

2FF 
Control Panel: OF 000 Program: (28.8) OF F 

26 2IJC 26 762F 

27K 
One or the other of these stops occurs at the end of each 
problem having Problem Parameter 8 = "+2+15^ n . A White 
Switch start will cause the routine to go on to the next 
problem. 

FF Stops 

I. Control Panel: FF 220 Program: (9*5) FF jkkF 

2.6 272 26 (7.004) 

297 
This stop signifies that the number i of polynomial 
coefficients read in or computed is impossible; that is, 
either i < or i > 50. Repositioning the data tape to 
the beginning of the next problem and White Switching will 
will cause the routine to go on to the next problem. 

II. Control Panel: FF 221 Program: (6k. 6) FF 5^5F 

26 2LK 2.6 (203) 

-23L 
This stop signifies that the number of extrema of 
E (x) found in [-1,1] is greater than 51, which is impossible, 
A White Switch start will cause the routine to go to the 
next problem 

III. Control Panel:: FF 222 Program: (211) FF 5^6f 

20 2JX 20 750F 

310 
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If the routine finds fewer than k+2 extrema of E (xj 
in [-1,1] during an iteration of Step II, it first prints 
out the coefficients of E^(x) and the extrema it has found 
if D > 0, and then comes to an FF 222 stop "because at 
least k+2 extrema are needed to determine a new error 
function E (x). The operator should now set the Black 
Switch to Obey and execute a White Switch start. A Black 
Switch start will now cause the routine to go on to the 
next problem, whereas a second White Switch start will cause 
the routine to replace N by N+10 and begin a new search 
for the extrema of E (x): in [-1,1]. If the number of 
extrema found using the increased value for N is now greater 
than or equal to k+2, the routine will Qbntinue on and finish 
the problem just as though no stop had occurred except 
that N has been replaced by N+10 for the remainder of 
the problem. If, however, the number of extrema found is 
still; .less than k+2, the routine will again stop on the 
FF222 order; and, just as before, the operator may cause 
the routine either to go on to the next problem or to re- 
place N+10 by N+20 and begin still another search for the 
extrema of E (x) in [-1,1]. 
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Section V EXAMPLES. 

The following seven examples illustrate how to prepare 

the data tape for and interpret the output from this routine. 

Under "Data Tape", Problems 1 through 7 are listed in order, 

just as they appear on the data tape. We note first that 

Problems 1 through 6 all have T = 'MO/fl^" and Problem 7 

has T = "+2+15V 1 ; therefore, if no failures occur and the 

Black Switch is set to Ignore, the routine will process all 

the problems one after the other, finally coming to an OF 

stop at the end of Problem 7« Looking first at the input 

data for Problem 1 under "Data Tape", we see that (l) the 

degree k of the preliminary and final approximating 

polynomials will be h, since m = "+h. n and 5 ="+1+152", 

(2) Parity ="+.", (3) N = "+4.", (k) intermediate results 

will be printed since D < 0, (5) the function which is to 

be approximated is the polynomial A (x) = x and (6) Option 

A-l is to be employed. Turning to the output data for 

Problem 1 under ."output", we first see Problem Parameters 1 

through 6 printed to 22 decimal places with a space after 

the first 11 places, and then the coefficients of the 

polynomial A. (x) = x in order. Next we see the coeffi- 

cients of the error polynomial E (x) = -.0703125+1. 05^875x - 

1.875x + x corresponding to the preliminary approximating 

polynomial P°(x) = ,0703125 - 1.05^675x 2 + 1.875* and a 

list of the seven extreme (x * E. (x ) of E. (x) in the 

r ' . 1 r 1 

interval [-1,1] , which are (-1, + .10938), (-.88508, - O 099722), 

..., (+1, + .10938). Next come the coefficients of E ± (x) and 

its extrema in [-1,1] followed by the final approximating 

polynomial, which is, to eleven significant figures, the 

polynomial f£(x) = +.0917327^3537 - 1.21^3267108x 2 + 

2.03086ll837x\ and by min |*£(x)| = +.0924059^2516, 

xe[-l,l] 

i a 

max |K(x)| = + .091732763537, and 5=k+l ' "T?l ' = + - 1095T5 ' 

X€[-l,l] 2 
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Looking now at Problem 2, we see that it is an Option 

B-l problem and therefore finds an approximation to the 

polynomial A (x) = x left in the machine by Problem 1. 

i 
This approximation will be of degree k as in Problem 1, 

but this time we have Parity = "+2.". The answers of 

Problems 1 and 2 are seen to agree to about 20 decimal 

places, which is about the accuracy to be expected with 

routine AJ» 

In Problem 3, we again approximate A^x) = x by a 
polynomial of degree k, using Option C-l, with results 
that are identical to those obtained in Problem 1 except 
that intermediate results are not printed. Note that the 
"+1+154" at the end of Problem 3 is not read in by the 
main routine but is read in by the interlude and stored 
at A7 address 1024 + 2(10) +2 = 10h6 for later retrieval 
by the main routine. 

Looking next at Problem k, we notice that m = "+." 
and 5 = "+.1095"; we will therefore approximate A ± (x) = x 
by a polynomial whose degree is such that an error no 
greater than .5 = .1095 is introduced in Step I (See Section 

III). The required degree turns out to te^, as we should 

10 a 
have expected from the fact that £ i r | _ + .109375 

r=2 ryC"!. 

and the results are again identical to those obtained in 
Problem 1. 

Problem 5 is an Option B-2 problem which first approxi- 
mates the polynomial A. (x) = x by a polynomial A (x) 
whose degree is such that an error no greater than 
5 = 0OOI96 introduced in Step I; that is, by the eighth 

degree polynomial A.(x) = .001953125 - .O9765625X + 

t 6 J 8 
.78l25x - 2.l875x + 2.5x . A.(x) is then approximated 

f 
by a polynomial P v (x) whose degree is such that an error 

no greater than 6 = .0976 is introduced in Step I. 
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Problem 6 is a Mode II, Option A problem . The two 
auxiliary subroutines calculate F(x) = x and F r (x) = lOx 
respectively whenever entered by the main iroutine; the 
results agree with those obtained in Problem 1 to about 
six decimal places. (Although one would not normally use 
Mode II to find an approximation to a polynomial of degree 
< lj-9, this has been done in Problems 6 and 7 for purposes 
of illustration). 

Next is a Dummy Problem which has no effect other than 
to change the End of Problem Constant at A7 address Ik from 
its normal value of "+1-1" to "+1-5". The End of Problem 
Constant will remain equal to this new value until it is 
changed again by means of another Dummy Problem. 

Finally, Problem 7 is a Mode II, Option B Problem. 
The auxiliary subroutines left in the machine are used to 
calculate F(x) = x and F f (x) = lOx . As can easily be 
seen, reducing the End of Problem Test Constant has caused 
the routine to make one more Step II iteration than in 
Problem 6, and thus to find a better approximation to x 
than in Problem 6. 



DAT E October 17, 19^1 
PROGRAMMED BY£ 
APPROVED BY 




nj 



+1. 
+ . 

+1+152 
+1*. 

-1. 

+ . 
+ . 
+ . 
+. 
+ . 
+ . 
+ • 
+ . 
+ . 
+ • 
+1. 
+1+15^ 

+2. 

+k. 

+2. 

+I+I52 

+k. 

-1. 

+1+153 

+1+15^ 

+3- 

+ . 

+1+152 

+h, 

+1. 

+2+153 

00866+ 

(1) 22(1) 50(1) 

26l6F 0+10F 

(2) 8+F 0-102*KF 
02(2) 8+1F 
0-1024F 0L2F 
88F Q-102^F 

8j(3) OOF 

(3) 226^9F OOF 
26866N 
+1+15 1 )- 



/ 
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DATA TAPE 

Problem Number (=l) 
m 

Parity 
.8 
N 
D 



' A.(x) 



T 



Problem Number 


(=2) 


m 

Parity 

6 

N 

D 

Select Option E 

T 


. 


Problem Number 


(=5) 


m 

Parity 

8 

N 

D 

Select Option C 


t 



Enter AT; b Q : = (-10, 0) 

Store at A? addresses 102*1- - 10*1-2 

Store 1 at A7 address lOMi- 

Read T from data tapej store T at 

A7 address 10h6 

Leave A7 

Return to main routine 






■Interlude 
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+4. 
+ . 
+ . 

+.1095 

+10. 

+1. 

+1+153 

+1+154 

+5. 

+. 

+ . 

+.00196 

+10. 

+1. 

+1+153 

+3+153 

+o0976 

+1+154 

+6, 

+4. 

+ . 

+. 

+4. 

-1. 

+5+153 
OO650+ 

(1) 22(1) 50(1) 
26l6F 8-698F 
0+10F &+1F 

(2) 87698F 03(2) 
8j(3) OOF 

(3) 22621F OOF 
00700+ 

(4) 22(4) 50(4) 
26l6F 8-698F 
0+9F 8+10F 

(5) 87698F 0.3(5) 
8j(6) OOF 

(6) 22501F OOF 
22 63ON 
+1+154 

+ . 
+ . 
+ . 
+ . 
+ . 
+. 

+5+153 
0094- 



1 



Problem Number (= 


V 


m 




Parity 
5 




N 




D 




Select Option B 




T 




Problem Number (= 


5) 


m 




Parity 
6 




N 




D 




Select Option B 




Select Option B-2 
8* 




T 




Problem Number (= 


6) 


m 




Parity 
8 




N 




D 




Select Mode II 





Enter A7; store (f ) at A7 address 698 

compute [ (f ) ] 

Leave A7 

Return to main routine 



Enter A7; store (f ) at A7 address 698 

compute 10 [(f)] 

Leave A7 

Return to main routine 



Auxiliary 
Subroutines 



T 



> Dummy Problem 

Changes A7 address l4 to "+1-5" 
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269N 
+1-5 

+7- 

+. 
+. 

+k. 

-l. 

+5+153 

OO65O+ 

2263ON 

+2+15^ 



11 



} 



Problem Number (= 7.) 

m 

Parity 

5 

N 

D 

Select Mode II 

Option JB 
T 
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OUTPUT 



Problem 1: 



+ 10000000000 00000000000 + 1 
+ 00000000000 00000000000 -155 
+ 1+0000000000 00000000000 + 1 



+ li-OOOOOOOOOO 00000000000 + 1 
+ 10000000000 ooooooooooo +152 
- 10000000000 ooooooooooo + 1 



+ ooooooooooo 
+■ ooooooooooo 
+ ooooooooooo 
+ ooooooooooo 
+ ooooooooooo 

+ 10000000000 

- 70312500000 
+ 105^6675000 

- 18750000000 
+ ooooooooooo 
+ ooooooooooo 

+ 10000000000 

- 10000 + 1 

- 88508 + 

- 53921+ + 
+ 00000 -155 
+ 5392*4- + 
+ 88508 + 

+ 10000 + 1 

- 91732763537 
+ 1211+3267108 

- 20308611837 
+ ooooooooooo 
+ ooooooooooo 

+ 10000000000 
- 10000 + 1 

- 89325 + 

- 55713 + 
+ 5^365 - 20 
+ 55713 +. 
+ 89325 + 

+ 10000 + 1 

+ 91732763537 

- 121^3267108 

+ 20308611837 
+ ooooooooooo 
+ ooooooooooo 
+ ooooooooooo 



ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 



-155 + ooooooooooo 

-155 + ooooooooooo 

-155 + ooooooooooo 

-155 + ooooooooooo 

-155 + ooooooooooo 
+ 1 



ooooooooooo - 1 

ooooooooooo + 1 

ooooooooooo ■+ .1 

ooooooooooo ^155 

ooooooooooo -155 

ooooooooooo + 1 



+ 10938 + 

J. 99722 - 

+ 79911 - 

- 70313 - 
+ 79911 - 

- 99722 - 
+ 10938 + 





1 
1 

1 
1 
1 





+ ooooooooooo 
+ ooooooooooo 
+ ooooooooooo 
+ ooooooooooo 
'■+ ooooooooooo 



ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 



ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 
ooooooooooo 



517281+11-6231 - 1 - 

0072003^572 + 1 + 

25685^5616 + 1 + 

ooooooooooo -155 + 

ooooooooooo -155 + 

ooooooooooo + 1 



+ 91733 - 1 

- 92356 - 1 
+ 921+06 - 1 

- 91733 - 1 
+ 92^06 - 1 

- 92356 - 1 
+ 91733 - 1 

51728^+6231 - 1 

0072003^572 + 1 

256851+65616 + 1 

ooooooooooo -155 

ooooooooooo -155 

ooooooooooo + 



13203^56937 
1+56071+30253 
ooooooooooo 
ooooooooooo 
ooooooooooo 



5711+2771+630 
72373388171 
ooooooooooo 
ooooooooooo 
ooooooooooo 



+ 132031+56937 
- 1+5607^30253 
+ ooooooooooo 
+ ooooooooooo 
+ ooooooooooo 



571^2771+630 
72373388171 
ooooooooooo 
ooooooooooo 
ooooooooooo 



■155 
-155 
■155 
-155 
-155 



-155 
■155 
-155 
-155 
-155 



- 19 
- 19 
-155 
-155 
-155 



- 19 
- 19 
-155 
-155 
-155 



+ 921+059^2516 8537^0300^1+ - 1 
+ 10937500000 ooooooooooo + 



+ 91732763537 517281+1+5977 - 1 



LOCATION 



781 



78^ 



626 



63O 



6^0 



6k6 



ORDER 



NOTES 



(210) 8J1(210) 83666F 
2^33F 22626F 

(211) FF5^6F 20730F 
OOF 22721F 
OO626K 

(0) 001F 50(0) 
26l6F 0K4F 

(0.01) 88F OS (jl) 

03(0.01) OK5E 
(0.02) 88F OS982F 

03(0.02) 88F 
8S2(J) 88F 
8S14F 852 (J ) 
804(J1) 82(0.1) 
(O.l) 8KF 8S2(J) 
8KF OK52F 

(1) OS(N) 03(1) 
9KLF 8L2800F 
8S(N) 8S52(N) 
1S4090F 1S^092F 
2K2^F 8j(l.l) 

(l.l) 22(1.1) 50(1.1) 
26l6F 8K2F 

(2) 8S(J) 81 (N) 
8S(N) 3KLOOOF 

(3) 93(3.1) 3KF 



From (9) 



From (26.1) 



From SADOI 



PAGE 3 



KA 2 



From 2(7) 



From 1(7) 



Leave A7, Go to 1(210 ); Go to 

(11) 

Stop if Obey, Then Bl. Sw. 

goes to (210), Wh. Sw. to (7.004) 

FF222J Wh. Sw. once, Stop if 

Obey 

Then Bl. Sw. goes to (27), Wh. 

Sw. to (26.1) 

BEGIN Program F 



Enter A7; b : = (-4, 0) 



Read k constants (+5+152, +50., 
+1., +5+153) 



Read 5 constants (+68-20, +2., 
+1.57076268,-. 6432292^.0727102); 
Read Chebyshev Print Constant. 



Read End of Problem Test Constant 



F: = Chebyshev Print Constant -1 



F: =0; b Q :=(-52, 0) 



F: =1, Refer to b ; 1^: 
b :=(-24, 0)| Leave A7 



= 2800 



2 i " 1 :=2 



(i)._ . + (i"2). 



t£""':= - t^~ ~ } % b:= (-1000,0) 
I b 3 :=(0,0) 
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Problem 2: 



+ 20000000000 00000000000 
+ 20000000000 00000000000 
+ 40000000000 00000000000 



+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 10000000000 

- 70312500000 
+ 10546875000 

- 18750000000 

+ 00000000000 
+ 00000000000 

+ 10000000000 



00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 

00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 



+ 00000 -155 - 70313 - 

+ 5392^ + +799H - 

+ 88508 + - 99722 - 

+ 10000 + 1 + 10938 + 



- 91732763537 
+ 12143267108 

- 20308611837 

+ 00000000000 
+ 00000000000 
+ 10000000000 



51728446062 
00720034572 
25685+65657 

00000000000 
00000000000 
00000000000 



+ 00000 -155 - 91733 

+ 55713 + + 92406 

+ 89325 + - 92356 

+ 10000 + 1 + 91733 



+ 91732763537 

- 12143267108 
+ 20308611837 

+ 00000000000 
+ 00000000000 
+ 00000000000 



51728++6062 
0072003+572 
25685+65657 

00000000000 
00000000000 
00000000000 



+ 92+059^2516 85374029790 
+ 10937500000 00000000000 



+ 1 
+ 1 
+ 1 

-155 
-155 
-155 
-155 
-155 

+ 1 

- 1 
+ 1 
+ 1 

-155 
-155 

+ 1 

1 
1 
1 



- 1 
+ 1 
+ 1 

-155 
-155 

+ 1 

1 
1 

1 . 
1 

- 1 
+ 1 
+ 1 

-155 
-155 

+ 

- 1 

+ 



+ 40000000000 00000000000 
+ 10000000000 00000000000 
- 10000000000 00000000000 



+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 



00000000000 
00000000000 
00000000000 
00000000000 
00000000000 



+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 



00000000000 
00000000000 
00000000000 
00000000000 
00000000000 



+ 1 
+152 
+ 1 

-155 
-155 
-155 
-155 
-155 



-155 
-155 
-155 
-155 
-155 



+ 91732763537 51728445469 - 1 



31 - 



Problem 3" 



+ 30000000000 00000000000 + 1 
+ 00000000000 00000000000 -155 
+ 40000000000 00000000000 + 1 



+ 40000000000 00000000000 + 1 

+ 10000000000 00000000000 +152 

+ 10000000000 00000000000 + 1 



+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 10000000000 

+ 91732763537 
- 12143267108 

+ 20308611837 

+ 00000000000 
+ 00000000000 
+ 00000000000 



00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 



-155 
-155 
-155 
-155 
-155 

+ 1 



51728446231 - 1 

00720034572 + 1 

25685465616 + 1 

00000000000 -155 

00000000000 -155 

00000000000 + 



+ 92405942516 85374030044 - 1 
+ 10937500000 00000000000 + 



+ 00000000000 00000000000 
+ 00000000000 00000000000 
+ 00000000000 00000000000 
+ 00000000000 00000000000 
+ 00000000000 00000000000 



+ 13203456937 
- 45607430253 

+ 00000000000 
+ 00000000000 
+ 00000000000 



57142774630 
72373388171 

00000000000 
00000000000 
00000000000 



-155 
-155 

-155 

-155 
-155 



- 19 
■ 19 
■155 
•155 
■155 



+ 91732763537 51728445977 - 1 



Problem 4; 



+ 40000000000 00000000000 
+ 00000000000 00000000000 
+ 10000000000 00000000000 



+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 10000000000 

+ 91732763537 
- 12143267108 

+ 20308611837 

+ 00000000000 
+ 00000000000 
+ 00000000000 



00000000000 
00000000000 
00000000000 
00000000000 
00000000000 
00000000000 

51728446231 
00720034572 
25685465616 

00000000000 
00000000000 
00000000000 



+ 92405942516 85374030044 
+ 10937500000 00000000000 



+ 1. + 00000000000 00000000000 -155 
-155 +1095000000004802132025 + 

+ 2 +10000000000 00000000000 + 1 



-155 
-155 
-155 
-155 
-155 
+ 1 

- 1 
+ 1 
+ 1 
-155 
-155 
+ 

- 1 

+ 



+ 00000000000 00000000000 -155 

155 
155 
155 
155 



00000000000 00000000000 
00000000000 00000000000 
00000000000 00000000000 
00000000000 00000000000 



+ 13203456937 571^2774630 - 19 

- 45607430253 72373388171 - 19 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 



+ 91732763537 51728445977 - 1 
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Problem 5* 



+ 50000000000 00000000000' + 1 ■ + 

+ 00000000000 00000000000 +155 + 

+ 10000000000 00000000000 + 2 + 



00000000000 00000000000 -155 

19599999995 989^9216306 - 2 
10000000000 00000000000 + 1 



00000000000 
00000000000 

00000000000 

00000000000 
00000000000 
10000000000 



00000000000 -155 

00000000000 -155 

00000000000 -155 

00000000000 -155 

00000000000 -155 

00000000000 + 1 



+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 
+ 00000000000 



00000000000 
00000000000 
00000000000 
00000000000 
00000000000 



-155 
-155 
-155 
-155 
-155 



+ 19551250000 00000000000 - 2 

+ 19531250000 00000000000 - 2 

- 97656250000 00000000000 - 1 
+ 78125000000 00000000000 + 

- 21875000000 00000000000 + 1 
+ 25000000000 00000000000 + 1 

- 17578125000 00000000000 - 1 
+ 5273^375000 00000000000 + 
-23^37500000 00000000000 + 1 
+ 28125000000 00000000000 + 1 
+ 00000000000 00000000000 + 

+ 19531250000 00000000032 - 1 

+ 19531250000 00000000000 - 1 



+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 



+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 



+ 195312^9999 99999999820 - 1 
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Problem 6: 



+ 60000000000 00000000000 + 1 

+ 00000000000: 00000000000 -155 

+ 40000000000 00000000000 + 1 

+ 3388131T890 172013562T3 - 20 

- 390908822^ 70*1-98371048 + 
+ 11723615957 22835053456 + 1 

- 10000 + 1 + 21855 +0 

- 84302 + - 13302 + 

- 4l04l + + 32718 - 1 

- 91752 - 20 - 33881 - 20 
+ 4io4i + + 32718 - 1 
+ 84302 + - 13302 + 

+ 10000 + 1 + 21855 + 

+ 73088277868 85538390285 - 1 

- 12156123013 89102779444 + 1 
+ 20694357456 51392011642 + 1 



- 10000 

- 89867 

- 55139 

+ 10765 
+ 55139 
+ 89867 

+ 10000 



1 



20 



1 



+ 73088 - 

- 97539 - 
+ 10781 + 

- 73088 - 
+ 10781 + 

- 97539 - 
+ 73088 - 



+ 92025754835 93762716152 
- 12139916368 34183746933 
+ 20299401271 62308492588 



- 10000 

- 89314 

- 55719 
+ 22353 
+ 55719 
+ 89314 
+ 10000 



1 



20 


1 



+ 92026 

- 92330 
+ 92098 

- 92026 
+ 92098 

- 92330 
+ 92026 



+ 92025754835 93762716152 - 
- 12139916368 34183746933 + 
+ 20299401271 62308492588 + 



+ 92330259242 66681706816 - 1 
+ 00000000000 00000000000 -155 



+ 40000000000 00000000000 + 1 
+ 00000000000 00000000000 -155 

+ 10000000000 00000000000 + 1 



- 71733102720 59855996485 
+ 31896084420 04474714304 



+ 39175273810 51139068191 
- 38381180422 46048411403 



+ 80755575252 03070932851 
- 27883017678 24373337840 



- 20 

- 20 



-20 
21 



20 
19 



+ 80755575252 02070932851 
- 27883017678 24373337840 



20 

19 



+ 92025754835 93762716152 - l 



- *- 



Dummy Problem : 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 00000000000 00000000000 -155 

+ 10000000000 00000000000 - 5 



+ 00000000000 00000000000 -155 
+ 00000000000 00000000000 -155 
+ 00000000000 00000000000 -155 



Problem 7» 



+ 70000000000 00000000000 + 1 
+ 00000000000 00000000000 -155 
+ 40000000000 00000000000 + 1 



+ 40000000000 00000000000 + 1 
+ 00000000000 00000000000 -155 
+ 10000000000 00000000000 + 1 



+^881517890 17201356273 -20-7173310272059855996^85 -20 

! 390SS 70^983710^8 + + 31896084420 044 T 47l4304 -20 

+ 11723615957 22835053456 + 1 

- 10000 + 1 + 21855 ■+ 

- 84302 + - 13302 + 

- 4io4i + + 32718 - 1 

- 91752 - 20 - 33881 - 20 

+ 4l04l ■+ + 32718 - 1 

+ 84302 + ' - 13302 + 

+ 10000 + 1 + 21855 + 



+ 73088277868 85538390285 - 
- 12156123013 89102779W + 

+ 206914-357*1-56 51392011642 + 



1 
1 

1 



+ 39175273810 51159068191 

- 38381180^22 46o484ll403 



- 10000 + 1 + 73088 - 1 



- 89867 

- 55139 
+ 10765 
+ 55159 
+ 89867 

+ 10000 



+ 
+ 
- 20 
+ 
+ 

+ 1 



- 97539 - 
+ 10781 + 

- 73088 - 
+ 10781 + 

- 97539 - 

+ 73088 - 



1 



1 



1 
1 



20 
21 



+ 92025754835 93762716152 - 1 + 80755575252 02070932851 
- S916368 34183746935 +.1 - 27883017678 24373357840 
+ 20299401271 62308492588 



20 
19 



+ 1 



- 10000 

- 89314 

- 55719 
+ 22353 
+ 55719 
+ 89314 

+ 10000 



+ 1 


+ 92026 - 


- 1 


+ 


- 92350 - 


■ 1 


+ 


+ 92098 - 


- 1 


- 20 


- 92026 • 


- 1 


+ 


+ 92098 - 


- 1 


+ 


- 92550 ■ 


- 1 


+ 1 


+ 92026 ■ 


- 1 



- 35 - 



+ 92161897888 6356081^550 
- 12l48o49122 23065368657 
+ 203048lll64 4579^152359 



- 10000 

- 89315 

- 55732 
+ 27619 
+ 55732 
+ 89315 
+ 10000 



1 


21 



1 



+ 92162 

- 92162 
+ 92162 

- 92162 
+ 92162 

- 92162 
+ 92162 



+ 
+ 

1 
1 
1 
1 
1 
1 
1 



+ 92161897888 63560814550 - 

- 12148049122 23065368657 + 
+ 203048lll64 45794152359 + 



+ 92161930797 35305528393 - 1 

+ 00000000000 00000000000 -155 



+ 99831329529 77386309035 
- 38501767846 13423017976 



+ 99831329529 77386309035 
- 38501767846 13423017976 



21 
20 



21 
20 



+ 92161897888 63560814550 - 1 
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Williams Memory: 

Addresses 0-2 
Address 3 
Addresses 4-8 
Addresses 9-13 



Addresses 
Addresses 
Addresses 
Addresses 
Addresses 
Addresses 
Addresses 
Addresses 
Addresses 
Addresses 



l4, 15. 
16-429 
430-749 
750-784 

785-791 
792-831 
832-839 
840-851 
852-863 
864-1023 



Temporary Storage 

S-parameter for Floating Accumulator 

Unused 

Routine to change End of Problem Test 
Constant 

End of Problem Test Constant 

A7 (excluding input-output) 

Program F, A, B, C, or D. 

Program E (Executive routine) 

Unused 

Yl 

Constants 

Temporary storage 

Problem Parameters 1-6 

Temporary storage and constants 



Drum: 



A7 Addresses 1024-1123 
(Drum addresses 2560-2659) 

A7 addresses 1124-1223 
(Drum addresses 2660-2759 ) 

A7 addresses 1224-1229, 
(Drum addresses 2760-2765) 

A7 addresses 1230-1425 
(Drum addresses 2766-2961) 

A7 addresses 1426-1529 
(Drum addresses 2962-3065) 



A7 addresses 1530-1633 
(Drum addresses 3066-3169) 

A7 addresses 1634-1735 
(Drum addresses 3170-3271) 

A7 addresses I736-I837 
(Drum addresses 3272-3373) 

A7 addresses 1838-1939 
(Drum addresses 3374-3475) 



Coefficients of A. (x) or A. (at) 
Coefficients of - (p£(x)) 
E.(-l), E.(0), E.(l) 
Unused 

Jli / \ 

fxl (r=l, • ••X) [the roots o f d ^ X in 

[-1,1]] if Parity = 0; orjJfeT (r=l .../£ ) 
if Parity = 2 l 

E*(x ) (r=l. ...t), the extrema of E,(x) in 
ir •*- 



n,i] 

Coefficients of E. (x) 



it'(x ); later, the quantities z, ..., z^ 
entering into the calculation of AE k (x) 
Unused 
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A7 addresses 19^0-204l 
(Drum addresses 3^76-3577) 

A7 addresses 20^2-21^5 
(Drum addresses 3578-3681 ) 

A7 addresses 2lk6-h095 
(Drum addresses 3682-5631) 

A7 addresses k096-kkl6 
(Drum addresses 5632-5952) 

A7 addresses kklJ-kjJ'T 
(Drum addresses 5953-6273) 

A7 addresses ^738-5058 
(Drum addresses 627^-659^) 

A7 addresses 5059-5379 
(Drum addresses 6595-6915) 

A7 addresses 538O-5576 
(Drum addresses 6916-7112) 

A7 addresses 5577-5773 
(Drum addresses 7113-7309) 

A7 addresses 689O-8189 
(Drum addresses 8^26-9725) 

Drum addresses 9726-10,999 

Drum addresses 11,000-12,799 



Coefficients of AEL (x) 

}x L (r=l, ...,c) Ethe roots of 



d^i (x) 



dx 



in [-1,1]] 
Unused 



Program A (including input-output routines 
from A7) 

Program B (including input-output routines 
from A7) 

Program C 



Program D 

A7 input -output routine 

First 196 words of Program C 

Chebyshev polynomial Coefficients 

Unused 
SAD0I 



LOCATION 


ORDER 


NOTES PAGE 4 KA 2 


6^7 


(3.1) 


OKF 9F11F 
8KF 8KF 
8N2(J) 8s4(j) 




b n :=(0 > 0); Set up Space 
after 11 digits 


650 


(3.2) 


82(3.2) 8Mf(Jl) 
8S2(J) 8F2F 
8F2F 85 (N) 


From 3(3.l) 


k(j): =-: (Chebyshev Print 

Constant | 

Chebyshev Print Constant := 

-4Ghebvshev Print Constant |+1 

F:=tW 




(3.3) 


33(5) 8K2.F 

32 GO 07(H) 

0052 (N) 0S52(N) 


From (6) 


J J-l J 




(5) 


93(5) 0750 (N) 
00 (N) os(n) 
86 (J ) 1S^09^F 


From l(3.3) 
From (3.3)^ 


J-l J 

_ . (i) . i-l 
Fs t: 7 -r 2 

3 






8KF lN^O^F 




660 




83(7) 85^ (J) 
82(5.1) 15^09^F 




F : — - | Chebyshev Print Cohsta 
F: =.t ( : i) ^2 i - 1 




(5.1) 


8922F 03(6) 




Print | n Q2 = n n + ^ 




(6) 


12(6) 92(3.3) 




n-, :=n-.+2 




(7) 


8K2F 87 (J ) 
8S(J) 32(3) 
22(2) 8j (7.002) 


From 2(5) 


2 1 " 1 : = 2 1 


66? 


(7*002 


) 26999F OOF 
26626N 
+5+152 
+50. 
+1. 
+5+153 

+68-20 




Go to (0) 
Constants 



LOCATION ORDER 

+2. 

+1.57076268 
-.6+32292 
+.0727102 

JL 

+1-1 

00626K 

626 (7.00+) 20(7.oo4) 50(7.00+) 
(7.01) 2616F ok6f 

8KF 8S992F 

(7.1) 88f os(ni) 

630 03(7.1) 8F2F 

8F2F 8F2F 

(7.2) 0K6F 05(N1) 
8922F 02(7.2) 
8F2F 8F2F 
J8F 8S10(J) 
8o(ji) 82(7.21) 

(7.21) 93(7.25) 8o(ji) 

8o(ji) 83(7.22) 
j8f 8sio(j) 

6+0 (7.212) 5KF 8KF 

(7.21+) 8S+(j) 55(N) 
8922F 5L2F 

85^(J) 8o(j) 
6++ 83(11) 8KLF . 



NOTES 



PAGE 5 



KA 2 



From 2(203), 
1(210), (9.5), 
3(2TA),(28) 



Coefficients of Rand approxima- 
tion to sine 



Chebyshev Print Constant 

End of Problem Test Constant 

BEGIN Program A 

BEGIN Problem. Stop if 
Obey. 

b Q : = (-6, 0) 

max. , : = 
l-l 

Read Problem Parameters 1-6 



From 2(23.4) 
From 4(7.21+) 



b Q :=(-6,0) 

Print Problem Parameters 1-6 



10(j): = First number of Problem 
Parameter 7. refer to b^ f 
F:=f-. 5x10152 . go to (7-21) if 
not Mode I, Option A I 
Go to (7.25) if Mode I, Option 
A; F^F-^xlO 1 ^ 2 J 

F:=F-..5xl0 152 ; go to (7.22) if 
Option C or Mode II | 
BEGIN Mode I, Option B. Read in 
,3xl0 1 53 or T 
T3 5 :(-5,0) 

4(j):=0 or 4(j)+l; F:= Polynomial 

coefficients already in machine 

Print: n c :=n 1 -+2 
5 5 

+(J): = +(j)-i 



L OCATION 



660 



670 
671 



ORDER 

8^(j) 93(7. 23A) 

(7.22) 8798^F 806(J1) 
83(8) 81(7.23) 

(7.23) 8511F ^-OF 
(7.2k) 26F 50(7.2^) 

26(A7) 0K1000F 

85102^ asio(j) 

(7.25) LK1023F 8S(J) 

Sl^^f 92(9) 

(8) 8J756F 03(8.1) 
88f 92(8.1) 

(8.1) 751026F 53(8.2) 

(8.2) asiO(j) 80(J1) 

(9) 83(210) 8EOJ 1 
Sk(j) 8S(J) 
802(J1) 82(9. *0 

(9.k) 93(10) 8J(9.5) 
(9.5) TF^kkF 26(7.00^) 

(10) 72(10) 8510(J) 
7S102^-F 7S(N) 
8922F 92(8) 

(11) 8KF 8S8(J) 
85(J) 8S2(J) 
8L2796F 6kf 
6f6if 8sMj) 
83(12) 92(9**0 

(12) 4KL000F ^-3(13) 



From 1(7.21) 



From (7.21) 



From 1(7.22), 
2(10) 



From (8) 



From 1(7.25) 



From ^(11) 



From (9*k) 



NOTES PAGE 6 KA~2~ 

153 
F: = 2F - 5x10 ^ J . 

Go to (8) if Mode II 

BEGIN Mode I, Option C. Read 
in Interlude 

b :=( -1000,0) 

10(j): = First polynomial 

Coefficient 

i:=-l, refer to t> 

iu : =^09^ = -2 mod ^096 ) 

Go to (329); Go to (8.1) if 
Mode I, Option C 



From 3(7.21*0 
(210) 



Read 1 Polynomial Coefficient; 
refer to Id,- ! 

10( J ) : =Polynomial Coefficient; 
F+=F -.5x101 52 J 

Go to (210) if finished reading 
Polynomial coefficients 
i: = i+1 

F:=F-50; Go to (9.^) if i > 50 



FF220; Reposition Data tape and 
Wh. S¥ a 
n : =n +2 



V 



102^+ru and (N) + n^ : =Polynomial 
Coefficient 
Print 

BEGIN Mode I, Step I. £ =k+1 l a r l : = 



k: = i 

n c : =2796 
5 

ng:=2i; lj-(j):=i 

Go to (9.k) if i < 

^:=( -1000,2) 



LOCATION 
672 



680 



69O 



(13): 



(1*0 
(15) 
(16) 



(16.05) 
(16.06) 

(16.07) 



(16.1) 



(18) 



698 



(19) 



(20) 



ORDER 

3K2F 8K1023F 
Qkk(j) 8SMJ) 
83(1^) 93(16) 
5F58F 32(15) 
93(13) ^2(13) 
852 (Nl) 802(J) 
82(23) 65(H) 
83(16.05) 61(B) 
57^09^ 8j(i6.o6) 
1937F L060F 
101F L^2S3 

if02S3 50(l6(,07) 
2616F 8^8(J) 
8s6(j) 856(Nl) 
8o6(j) 82(16.1) 
92(23) 856(J) 
8S8(J) 33(18) 
3K2F EL2F 
8KF 6n(n) 
82(21) 2KF 
2F60F 852(j) 
QBk(j) KI^F 
JL2F 8K1022F 

ekk{j) 8SMJ) 
83(20) 8kf 
6s(H) 93(22) 
6l(N) 57^09^F 



NOTES 
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From (15) 



From 2(13), 
(23) 



From l(l6.l) 



From (21) 



From 3(19) 



■b 3 :«(-2,0) 
MJ):=MJ)-1 



^y=^\ 



F:-=?m-]& 

Go to (23) if finished with 
Step I; F: = a|. k ) 

F: = K ; | • t k 

A: = 2~ 3 - n 6 

A: = A/2 + (Exponent of F) 

(Exponent of F): -A 



V -|^)| • t*/ 2 ^i +1 |aW| 

F: = S - F^. 

Go to (23) if finished with Step I 

b used as binary switch 

^3: =(-2,0); n^:=n^_ 2 

f:=-|^ k) | 

If \ =0, bypass 

BEGIN Subtraction of T,(x) • r^: =n^ 



k(j): = k; n 2 :=n 2 -^ 

n 5 :=n 5-2 

MJ): = k(j) -2 






LOCATION 



ORDER 



699 
700 



710 



720 



725 



2k (n) 2S(N) 

(21) 92(19) 5J58F 

(22) FI2F 8KL023F 
Qk2 (J) 8S2(J) 

(23) 93(16) 8K7F 
87 (Jl) 8010 (J) 
83(23.1) 8J(205) 

(23.1) 852 (J) 8S(J) 
03(23.2) 88F 
8s6(Nl) 88F 
8S10(J) 93(23.3) 

(23.2) 751028F 8S6(N1) 
751030F 8S10(J) 

(23.3) 8F2F 858(J) 
8922F 8F2F 
7KF 8KF 

(23.^) 8S4(J) 75 (N) 

7S1024F 85^ (J ) 
8o(j) 83(7.212) 
8KLF 8^(J) 
7I2F 93(23.^) 

(26) L5(27-l) ^22(204) 

(26.1) 26(211) 50(26.1) 

26l6F 8KL0F 
8^8(Nl) 8s8(Nl) 
5KF 8KLF 

(26.2) 8S4(J) 55163^ 



NOTES 
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From 1(18) 
From Ml9) 



From l(l6), 
(16.1) 



From 1(23.1) 



From 3(23.1) 



From kfej.k) 



From 2(20^+) 
From 1(211) 



From (27) 



r r k r 



ry^tyn^ 



n 6 s=n 6-2 
k: = k-1 



^53 



153 



F:=35xHT -(.3 x 10 w or T) 

Go to (23.1) if Option A-2, 

B-2, or C-2 

BEGIN Modified Option o i:=k 

Binary Switch; Go to (23.2) 
if Option B 
o:=S ; Read T 



10(J): = T 

f.8 =Z i .^[a^ 1 " 1 

r=j+l l r ' 
Print 



M«l): = or k(j) + 1 



Go to (7.212) for second itera- 
tion of Step I 



BEGIN New Step II Iteration with 
N:=N+10. Go to (211 ) for FF222 
Stop) 



N: = N+10 



b^: = (0,0) 

4(j): = 1 or k(j) + 1 



L OCATION 
726 



730 



7^0 



7^9 



750 



ORDER 

5S(N) 5L2F 

85976'F 80^ (J ) 

83(26.3) 8J(206) 
(26.3) 8KLF 8l4-i+(j) 
(27) 93(26.2) 50(27) From (211) 

26l6F 92(27.35) 

(27.1) 001F 50(27.1) From 2(20*0 
26l6F 85 (J ) 
8S4(J) 8F2F 
8F2F 2KF 

(27.2) 25102^F 2016 3^F From 1(27.3) 
8922F 23(27^3) 

(27.3) 85^(J) 80MJ1) 
8sMj) 83(27.2) 
851022F 8F2F 
8922F 851018F 
8922F 858 (J) 

(27.35) 8922F 2K51F From 1(27) 
(27 .h) 25102^F 2S(N) 

23(27A) 8K1021F 

876 (Jl) 81«L0 (J) 

83(27.5) 8J:{7.00^) 
(27.5) 8J(28) OOF 
(28) OFF 26 (7.004) 

2675ON 

00750K 

J0^30F 50750F 



NOTES 
n 5 :=n 5+2 
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Go to (206) for new Iteration 
of Step II 



BEGIN Printout of Results 



^(j): = i(or j if Modified Option) 



F:^ 1 -e f (Final Approximating 
r Polynomial) 



Print; n 2 :=n 2 +2 



MJ):=MJ)-1 

F:=max \W. (x) 

xe[-l,l] 
F:=min 



or |E. (x) 
1 J 



^i(or J) |a .(r) /2 r-l 



r=k+l ' r 
BEGIN Going on to next problem 



FJ=-.15 x 10 15 +T 

Go to "begin next problem if 
T=.lxl0 1 ^ 



OF Stop; Wh. Sw. to begin next 

problem 

Go to (200) 



Record A7 Input -output routines 
on Drum at 6916-7112 J 



LOCATION 



ORDER 



751 26 (yi) 006916F 

752 OOI96F 26999F 
2675ON 
OO626K 

626 (28 ,5) OOIF 50(28.5) 
26l6F 8KF 

8n8(j) 82(28.6) 

(28.6) 93(28.9) 8KL021F 
630 876 (Jl) 8410 (J ) 

85(28.7) 8J762F 

(28.7) 8j(28.8) OOF 

(28.8) OFF. 26762F 

(28.9) 2KF 2KF 
(29) 85 (J) 8SMJ) 

21 (N) 2S1124F 
241024F 2S(N) 
(29.1) 22(29.1) 85MJ) 
80^ (Jl) 82(29) 
6ko 8KF 8s6(ni) 

8nMni) 83(30) 

8kk (Jl) 8S1004F 
5KF 3K1000F 
83(29.2) 85 (J) 
8S976F 852 (J ) 
8S9T8F 92(29.3) 
Gkj (29.2) 5L2F 85 (J ) 



NOTES 
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From 2(205) 



From (28.6) 
From 1(29.1) 



BEGIN Program B 



From 6(29.1) 



FMZ 1 la^^ll 
l ¥=k+l ' r ' ' ' 

Go to (28.9) if m < i 
F:=~.15xl0 15 + T 



Go to (2033 to begin next problem 
if T=.lxl0 154 



OF Stop. Wh. Sw. to go to (205), 
next problem I 

BEGIN First Iteration of Step II 



MJ): = i 

122kF + n : = - P^ 
2 r 

(N) + n 2 : = .e* 

n 2 ;=n 2+ 2 

Go to (29) if all coefficients 

transferred 

5: = 

F: = -[Parity |f Go to (30) if 
Parity = 
BEGIN Deleting Coefficients if 
Parity f 



Go to (29.2) if Parity 
976F: = i 
978F: = k 



= 1 



LOCATION 



648 



65O 



660 



670 



ORDER 



NOTES 
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67k 



8o4(jl) 8S976F 
852 (J ) 8o4(«Jl) 

(29.3) 8S978F 85976F 
86984F 8S976F 
85978F 86984F 
8j(29.4) 93(29.8) 

(29.4) 15(29.7) L02S3 
42(29.5) 42(29.6) 

(29.5) L5S3 10F 

(29.6) :.50402F OOF 

40S3 2633F 
(29 o 7) OOF 00550F 
(29.8) 8S978F 81976F 
(29-9) 8s4(j) 55 (N) 
3S(n) 5L4f 
854 (J) 83(30.05) 
844 (ji) 33(29.9) 
(30) 85 (J ) 8S976F 
852 (J ) 8S978F 
3KF 3F61F 

(30.05) 8510 (Hi) 82(30.12) 

2KF 85976F 
8S4(J) 8F2F 
8KF 8n4(N1) 
82(30,06) 8F1F 

(30.06) 93(30.07) 8F2F 

(30.07) 25 (N) 8922F 



From 8(29.1) 



From 1(29.6) 



By 1(29.4) 
By 1(29.4) 



From 5(29.3) 
From 3(29.9) 



From 3(29.1) 



From 2(29.9), 
1(207) 



From (30.08) 



976: = i-1 

978: = k-1 
i: = i/2 
Fs = k/2 

A: = 550 - Exponent of k/2 



Q: = 

Go to 3(29.3) 

k: = [k/2] 

4(j): = -i or 4(j) + 1 

Delete unneeded coefficients 

Go to (30.05) if finished 
deleting coefficients 

976F: = i 

978F: = k 

ly = (0, n ? ) 

BEGIN printout of^polynomial 
coefficients of P?(x). Go to 
(30.12) if D > 

4(j): = i 

Go to (30.06) if Parity = 



F: = E ; print 



LOCATION 



675 



680 



69O 



700 
701 



ORDER 



NOTES 
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2L2F 854 (J) 
8o4(jl) 8s4(j) 

(30.08) 83(30.07) 22(30.09) 

(30.09) 92(30.12) 8J(30.10) 

(30.10) L5(30.1l) 421(207) 
1934F L4755F 
42755F 22(42.62) 

(30.11) OOF 00(30.05) 

(30.12) OOF 8KF 

8n4(ni) 83(30.13) 

844 (Jl) 82(30.13) 

(30.13) 85(N) 8S1226F 
SfCLF 8S1002F 
2K1F 2K1F 
859T6F 8s4(j) 
5KF 5F57F 
55(n) 8sioo4f 
854(j) 8o4(ji) 
83(33) 93(34) 

8s4(j) JL2F 

831004-F 871002F 
54(n) 92(32) 
851004F 23(38) 
8S1228F 8KF 
8n4(ni) 83(36) 
844(ji) 83(35) 
851004F 93(38) 



(3D 



(32) 



(33) 



(34) 



From (30.05) 



From 1(30.12) 



From 1(36) 



From 2(33) 



From 2(32) 



4(J): = 4(J) - 1 



Go to (3^.09) Right if printing 
because C < k+2 j 

Go to (30.12) if printing he-j 
cause D < ! 

A: = (30.05) { 



2 (I28);=l6(71. 5l)l go to (42.62) 
to print extrema 



Go to (30.13) Left if Parity = 

Go to (30.13) Left if Parity = 2 

1226F: = E. (0) 

1002F: = 1 

Set binary switch = - 1 

4(J): = i 



t> 5 :=(0, n 3 ) 



V 



BEGIN Evaluation of E. (1002F) 



Go to (34) if finished evaluating 



4(J): = 4(J) = 1; n 5 :=n 5 -2 



0/ 
F:= E? (x); go to (38) if x=-l 



1228F: = E. (1) 

Go to (36) if Parity = 

Go to 4(34) if Parity = 2 



LOCATION 



702 



710 



7L^- 



^30 



kko 



ORDER 



(35) 8llOO^F 93(38) 

(36) 811002F 8S1002F 
2K2F 93(31) 

(38) 8S122^F 93(^2.62) 

(te.62) 81(208) 50(^2.62) 

26l6F 8F2F 
5KF\ 8K1023F 
8K1023F 8^-980F 
(^2.63) 8Bif(j) 5520^2F 
895F 551530F 
895? 8K1023F 
Shk(j) %2F 
83(^2.63) 8J753F 

26753N 

00^30K 
(1+2.65) OOIF 50(^2.65) 
26l6F 5K1000F 
5K1000F 8K1023F 

(I4-2.8) 8bMj) 80976F 
83(^2.9) 55(N) 
5Sl63^F 8kif 
Qkk(j) 53(^2.8) 

(^2.9) 5KF 8KF 

8sM«J) 85(N) 

(1+3) 8b6(j) 80.F 



NOTES 
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From 3(3*0 
From 2(3*0 



1002F: = -1 



Set binary switch = -2 



y 



From (3k),k(3k), 122^F: = E. (-1) 

fe) - - , 0/ 

Go to (208) to continue iteration 

of Step II 

BEGIN printout of extrema 

\. = (0,0) 



From 2(30.10) 
2(127) 



From ^(^2.63) 



From 2(208) 



From 3(1*2.8) 



From 1(^-2.8) 



From 3(^3.3) 



b(j): = t-1 or k(j) -1; F: = x 

Print: F: = E. (x J 
y 1 x r 7 

Print 



n,_: = n n +2 
5 5 

Go to (128) if finished printing 

extrema 

Go to (200.5) 

BEGIN Program C 



BEGIN transfer coefficients of E 
from (N) to JL63kY, where E(x) = 

E*(x) or P*(x) 

1 v ' k 

\{J): = -1 or Ij-(j) +1 

Go to (^2.9) if finished 

l63^F + n c : = e 
p r 

n,-: = n._+2 
5 5 

BEGIN compute coefficients of E 1 
k(j): = 



MJ): MJ) + 1 



LOCATION 

kki 



ORDER 



k^O 



I46O 



From (^3.1) 
From 1(^3.1) 



k6l 



8KF 8N^(Nl) 

82(^3.1) 8^(J1) 
(1*3.1) 83(^3.2) 85MJ) 

2KL000F 93(^3.3) 
(ltf.2) 2KF 85MJ) 
(11-3.3) SL2F 571636F 

3S(N) 85MJ) 

80976F 5L2F 

82(li-3.5) 92(^3) 
(43.5) 3F59F 3F59F 
&(kk) 93(^5) 

(kk) L559F 0019F 

14-6(50) 2633F 
(I45) 8K1023F 8S1000F 

8512214F 8S1002F 

8KF 8n6(ni) 

83(145.06) 851000F 

8j(l45.02) 001F 
(1^5.02) 001F 50(145.02) 
(I1-5.0I4) 26(71.18) 50(145. Ok) From l(Tl.l8) 

26l6F 801002F 

8S1002F 8KF 
(I45.06) 8KF 8S980F From 3(^5) 

I4KF 8Nli-(Nl) 

83(^5.25) 8KF 

8S1000F 8S1008F 

851226F 8S1002F 



From l( i 43«5) 
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Go to (143.1) Right if Parity = 
Go to (14-3.1) Right if Parity = 2 



ru: = iu-2 

(W) + r-1: = r e 



Go to (^3.5) if finished 

BEGIN search for extrema in 
[-1,1] • n 3 : = i 



Go to, 1(1+3.5) 



Vr * L1 

E(-l): = E*(-l) 



Go to (^5.06) if Mode I 
F: = -1 



Go to (71.18) to compute F(-l) 



E(-l): = F(-l) - E(-l) 
t: =0 



Go to (145.25) if Parity = 

x _ : = 0: s : = 

n-1 ' 

E(0) 



LOCATION 


ORDER 
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14-68 


851636F 22(45.1) 
(14-5.1) 856 (J) 8S1010F 




Go to (45.1) right if Parity 


= 2 


470 


8K1F 92(45.25) 










(^5.25) 818 (Ml) 8S99^F 


From 2(45.06) 


h:=-N if Parity = 0, 1 if Parity 
= 2 




(45 o5) 868 (Hi) 8S996F 


From 5(69) 


F: = h/N 






87996F 8S998F 










8799OF 8I4988F 










87998F 84986F 










(45 07) 87996F 8S1020F 


From (69.I) 


s(h)-:- = sin(it/2 h/N) 






8KF 8n4(N1) 










83(45.75) 851020F 




Go to (45.75) if Parity = 






871020F 8S1020F 






s: = s*" 




48o 


fr5.75) 8J(46) 92(47.5) 


From 2(45.7), 
(*7) 








(46) L5(47) 42(51.8) 


Exit at (51.8): = (52) 






(47) 2633F 00(52) 




Go to (45.75) 






(47.5) OOF 851020F 


From (45. 75) 


Fs ■= s(h) 






(48) 8S1012F 8KF 


From 5(52.5) 


BEGIN evaluate E f (x)$ x: = F 






(50) 5KF 871012F 


By 1(44) 


*y = (-i, 0) 






5^(N) 52(50) 










8S1014F 8KF 




1014F: = E' (x) 






8n4(N1) 82(51 = 5) 




Go to (51.5) Right if Parity = C 




23(51-5) 851014F 




Go to (51.5) Left if Parity = 


= 2 


490 


871012F 87984F 
8S1006F 8j(50.l) 










93(60.5) OOF 
(50.1) L5(50.2) 42(63) 


From (50.2) 
From 6(50) 


Go to evaluate E(x) if Parity 
or 2 


r + 


k9k 


(50.2) 2633F 00(50.3) 




Go to 7(50) 





LOCATION 



1*95 



500 



510 



520 

521 



ORDER 



(50.3) 841006F 8S101^F 
92(51.5) OOF 

(51.5) 2KL000F 8KF 

8n6(ni) 82(51.8) 
851012F 8j(51.6) 

(51.6) OOIF 50(51.6) 

(51 . 7) 26700F 50(51-7) 
2616F 801014F 

(51.8) SSIOI^F 93F 

(52) 851014F 8S1022F 
818 (Nl) 8099^ 
83(68) 8KF 

8nioi4f 82(58) 

851020F 8S1016F 
851022F 8S1018F 
871010F 83(68) 
(52.5) 651010F 801018F 
8S1012F 851008F 
87IOI8F 8S101*4-F 

85IOI6F 87IOIOF 

801^F 861012F 
8J(53) 93(^8) 

(53) L5(5*0 ^2(51.8) 
(5*0 2633F 00(55) 
(55) 8KF 8N1014F 

82(58) 85IOO8F 
801016F 8S^(J) 



From (63) 



From 3(50)5 
M50) 



From Auxiliary- 
Subroutine 



By (46), (53) 
From 1(51.5) 



From (58 ) 



From (5^) 



From (51.8) 
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b : = (-1000,0) if Parity / 
Go to (51.8) if Mode I 



F: = x 



Go to Auxiliary Subroutine to 
compute F' (x) 



E' (x): = F r (x) - E' (x)j EXIT 

F: =-N-h 

Go to (68) if h = -N 

Go to (58) if E'(s) = 



Go to (68) if E f (s(h)). E f (s(h-l)) 
> 

BEGIN search for Root in 
[s(h-l), s(h)] 



F! = Z n + l= (z 2)n E 'K)n-( Z l)n ET ( Z 2) 



n 



E' (z ) - E' (z ) 
Go to (48) to evaluate E'^z^) 

Exit at (51.8): = (55) 

Go to 5(52.5) 

Go to (58) if E'(z n+1 ) = 
MJ): - (z l)n - ( z 2 ) n 



LOCATION 



522 



530 



540 



548 



ORDER 
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85982F 8N4(j) 
82(58) 851014F 
871010F 83(56) 
851010F 86984F 
8S1010F 93(57) 

(56) 85IOI6F 8S1008F 
851018F 8S1010F 

(57) 851012F 8S1016F 
851014F 8S1018F 

(58) 93(52o5) 93(58.5) 
(58.5) 8J(59) 93(60.5) 

(59) L5(6o) 42(63) 

(60) 2633F oo(64) 
(60.5) 85976F 8s4(j) 

5KF 5E57F 

(61) 5516 34-F 8S1014F 
854 (J) 8o4(jl) 
83(62) 93(62.5) 

(62) 8s4(j) JL2F 
851014F 871012F 
5^16 3^F 92(6l) 

(62.5) 8KF 8N6(Nl) 
83(63) 851012F 
8J(62.6) 001F 

(62.6) 001F 50(62.6) 

(62.7) 26(71.18) 5P(62.7) 
26l6F 801014F 



From 5(55) 



From. JJ (55) 



From 3(52), 

1(55), 4(55) 
From (60), 
4(70) 



From 7(50), 
(58.5) 



From 2(62) 



From 2(6l) 



From 1(71.18) 



Goto (58J if \(* 1 ) n -(z 2 ) n \ 

Go to (56) if E' (z -,) . E' (z ) 
> 



E, K)n + l :=E, K)n^ 2 
E, < Z A+1 : =E,(z 2 } n 



^ Z 2^n+1 : " Z n+1 



J 2 y n+l 



n+1' 



Go to (52.5) to iterate search j 
EXIT 



Exit at (63): = 64 

Go to (58.5) 

BEGIN Evaluate E(x). 4(j): 



= 1 



n 5 : = ^°' n 3^ 
1014F: = E(x) 



Go to (62.5) if finished 
4(J); = 4(J)-1 



Go to (63) if Mode I;=F: = x 



Go to (71.18) to compute F(x) 



LOCATION 



5^9 
550 



560 



570 



575 



ORDER 



8S1014F 851014F 

(63) 851014F 93F 

(64) 871002F 82(65) 

(64.1 ) 5KF 851000F 
4S2042F 4S2042F 
4S1426F 8KF 
8n4(N1) 82(64.4) 
8KF 8N1000F 
8S1018F 82(64.3) 
8K1F 841000F 

(64.2) 86984F 8sioi6f 
85100OF 861016F 
84ioi6f 86984F 
8S1018F 851018F 
801016F 82(64.3) 
851018F 92(64.2) 

(64.3) ■ OOF 851018F 

4S2042F 854 (Nl) 
8N984F 82(64.4) 
851018F 871002F 

(64.4) 93(64.5) 851002F 

(64.5) 4S1530F 8KLF 
8498OF 8S98OF 
8K51F 8098OF 
83(64.7) 8J(64.6) 

(64.6) FF545F 26(203) 

(64.7) 4L2F 53(71*02) 



NOTES 
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By (50.1), (59) 
From 1(62.5) 



From (71) 



From 5(64.2) 



From 5(64.1), 
4(64.2) 



From 3(64.1), 
2(64.3) 



From 2(64.5) 



E(x): = F(x) - E(x) 
F: = E(x); EXIT 
F: = E(x). E(x r ) 

F: = x 



Go to (64.4) if Parity = 



BEGIN calculate 



V-l 



1018F: = -|x J; go to (64.3) 

if x 1 = r " 

y :=r 



1018F: = y, 



j-^y.-y. x ; Go to (64.3) if 



- 3- J 
finished 



2042F+n^s=y J (=^T^) 

Go to (64.4) if Parity = 2 

F: = fx T . E(x , ) 

1 r-1 v r-1' 

F: = E(x r-1 ) 

1530F + n^: = E(x r _ 1 ) 

C : = L +i 

Go to (64.7) if ^ > 51 

FF221. Wh. Sw. to begin next 
problem 

n^:=n^+2; go to (71.02) if storing 
last extremum j 



LOCATION 


ORDER 
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576 


(65) 


93(67) 851002F 
83(66) 811002F 


From (6k) 






(66) 


8N101^F 83(68) 




F:=lE(x r _ 1 )| - |E(x r )| 




(67) 


851012F 8S1000F 


From (65) 


X r-1 : = X r 


58O 




85101^F 8S1002F 




E(x r _ 1 ): = E(x r ) 




(68) 
(69) 


85994F 808 (Nl) 
82(70) 851020F 
8SIOO8F 851022F 
8S1010F 8KLF 
&k-99kY 8S99^F 
808 (Nl) 83(69.1) 
8599^F 930*5.5) 


From 2(52), 
6(52), (66) 


F: = h-N 

Go to (70) if finished finding 

extrema 

s(h-l): = s(h) 

E(s(h-l)): = E(s(h)) 

h: = h+1 
F: = h-kL-N 




(69.1) 


8KLF 92(^5.7) 


From k(69) 






(70) 


OOF 80^(J1) 


From (69) 


F: = h-N-1 


590 




83(71) 8KLF 
8^99^-F 8S99^F 




t : -4+1 






8KLF 8KLF 
8S1012F 93(58.5) 








(71) 


5K2F 92(6^.1) 


From 1(70) 


Go to (6*1-. l) to store last 
extremum 




(71.02 


) 8J(209) OOF 


From (6k cj) 


Go to Program D 




(71.06) 3KF 6KF 


From 3(71- i 


BEGIN continue to set up Mode 






3f6if 6f6if 




Ta : = (0, n ? )| h 6 : = (0, n ? ) 






8KLF 8te(Nl) 










8S101^F 8698^ 




lOl^F: = m+1 


600 




8S1016F 2KF 
812 (Nl) 861014F 




1016F: = (m+l)/2) 


602 




8S101^F 8KF 




1014F: = -m/m+1 

.. — 



| LOCATION 


ORDER 
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603 


(71.08) 8S1018F 861016F 


From (71.1^) 


t: = or t+1 




84l01^F 8S1020F 




x: = -1 + (2t+l)/(m+l) 




871020F 851022F 








8799OF 8^988F 








871022F 8^986F 








871020F 2S1^26F 




1426F + n : = sin it/2 x fc 




2S2042F 8J (71.10) 




2042F + n : = sin it/ 2 x fc 


6lO 


(71.10) OOIF 50(71.10) 








(71.12) 26(71.18) 50(71.12) 
26l6F 2S1530F 


From 1(71.18) 


Go to (71c 18) to compute F 

;(sin it/2 x, ) 

1530F + n : = F(sin */2 x k ) 




851018F 802 (Hi) 








82(71.110 2L2F 




Go to (71.1*0 if finished 




8K1F 841018F 








(71.1*0 93(71.08) kKF 


From 3(71.12) 






4F56F 8kf 


■ 


n^: = n 2 




(71.16) 2K50F 2S102^F 








22(71.16) 8J759F 




Go to (130) 


620 
621 


(71.18) K5F ^2(71.20) 
(71.20) 2665OF 22F 
00626K 


From (^5-04) , 
(62.7), (71.12) 
From Auxiliary 
Subroutine 


BEGIN Go to Auxiliary Subroutn 
for F(x) 

BEGIN Set up Mode II 


626 


(71.3) 50430F 50(71.3) 

26 (yi) 006916F 

00196F OOIF 


From 2(129) 


Playback A7 Input-Output Rout: 




8511F 40F 




Read in Auxiliary Subroutines 


63O 


(71.32) 26F 50(71.32) 


From SADOI 






26 (AJ) 88F 




Read T 


632 


8S10(J) 8J(71.33) 




_ — _ . — 



LOCATION 
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633 


(71.33) 50^30F 50(71.33) 
26 (Yl) 007113F 
OOI96F 001F 

(71.34) J0430F 50(71.3^) 


a 


I Playback Program C ^30-626 




26 (Yl) 0062 7^F 


1 Record new Program G on Drum 




(71.36) 00320F 50(71-36) 








26 (A7) 8K2F 






640 


876 (Jl) 8S6(N1) 
0K51F 8KF 
(71.38) OS1634F 03(71.38) 




5: = 2T 




852 (Hi) 8S(J) 




i: = m 




8S2(J) 8S976F 




k: .=fflj i*. = m 




7KF 8S978F 




b : = (0,0); k: = m 




(71.4) 804(J1) 7L2F 




F: = m-1; n 7 : = 2m 




80^ (Jl) 82(71;^) 




s 




8KF 8S8(J) 




Z\ja (r) |: =0 

r=k+l' r ' 


6k9 


93(71.06) OOF 




Go to Continue setting up Mod« 




26756N 




Go to (201) 




00750K 






750 


J0430F 50750F 


Fronv-SADDI 


Record Program C 43O-626 on 




26 (Yl) 007113F 




Drum at 7H3-7309 


752 


OOI96F 26999F 
2675ON 








00430K 




BEGIN Program D 


^30 


(71.5) 001F 50(71.5) 


From 2(209) 




kjl 


26l6F 8510(11) 







LOCATION 



k32 



kko 



k^O 



458 



ORDER 
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8510 (EL) 82(71.52) 

(71.51) 8J750F 50(71.51) 

(71.52) 26l6F 8K1022F 
8il-980F 80978F 
82(71.57) 8510 (Kl) 
83(71.53) 8J(71.56) 

(71.53) 2K2F 2L768F 
85976F 85976F 

(71.5*0 8j(71.55) 002(30.05) 
(71.55) L5 (71.5*0 ^21(207) 

26(206) OOF 

OOF OOF 

OOF OOF 

OOF. OOF 

OOF OOF 

OOF 00(26) 

OOF 00(71.51) 

OOF L515(71-5l)' 

^2755F OOIF 



From 2(128) 
From 2(71.5) 



From 2(128) 



Goto (71.52) if D > 

Go to (127) to Print Extrema 

F: = 2-k-2 

Go to (71.57) if ^ -k-2 > Q 
Go to (71.53) if D < 
b 2 : = (0, 768) 

F: = i 

1(207): = 2(30.05) 



=Go to print Polynomial Coeffi- 
cients, Extrema 



Restore 2(^28); = (71-51) 



(71 ."56) L51^ (71.51) ^22(20*0 From 3(71*52) 2(20*0: = (26) 



(71.57) 26(20*0 NL2F, 
8K10F 8S1020F 
8KLF 861020F 
8S1020F 8KF 
8S1022F 856 (Jl) 
8S1018F 2KF 
8K1023F 84980F 



From 2(71,52) 



Go to (20*0 to playback Program 
A, FF222: Stop | 

BEGIN search for Max* and Min. 
Extrema 



1020F: = .01 
I max I : - 
1018F: = .1 x 10 155 



LOCATION 
459 

46o 



470 



480 



From l(71»6o) 



From 1(71.64) 



485 



ORDER 

(71.60) 8s4(j) 851022F JFrom 1 (71-70 ) 

2N1530F 83(71.64) 

251530F 82(71.62) 
(71,62) 211530F 8S1022F 
(71.64) 851018F 2N1530F 

83(71.66) 93(71-70) 
(71.66) 251530F 82(71.68) 
(71.68) 211530F 8S1018F 
(71o70) 2L2F 854 (J ) 

804(J1) 83(71.60) 

851022F 801018F 

8S1016F 8KF 

8N1016F 861022F 

8S1016F 8KF 

85i4f 84ioi6f 
82(71.7^) 92(72.6) 

(71.74) OOF 8j(71.76) 
(71.76) 24(71.78) 22(71.-78) 
(71.78) 22(71.82) 50(71.78) 

2616F 92(72.6) 
(71.82) OOF 50(71.82) 

2616F 8kf 
8kf 8n4(ni) 
83(72.5) 8KL023F 
84980F 8S4(J) 
2KF 2F58F 
2F58F 5KF 
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/^ 



4(J): -^-1 or 4(J) -1 
F: = |max| - |extremum| 

I max I : = |extremum| 
F: = iminl - lextremum 



Imin : = extremum; 



l6F: = I Imax! - (mini 



1016F: s- I Imax I - |mln| I7 |max 



F:=End of Problem Constant 

-m - in -hi 



Stop if Obeyj Bl. Sw. go to 
(71.78) left 
Wh. .Sw. go to (71.78) right 9 then 
(72.6) ' 



BEGIN Problem Termination Routine 



Go to (72.5) if Parity = 
4(J): =^-1 



V = (0, i\) 



LOCATION 


ORDER 
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14-86 


5F58F 8KLF 

804 (Nl) 8S1004F 




ty = (0, n 5 ) 




8S1004F 8KL022F 




1004F: = Parity -1 




844(Nl) 82(71.86) 




Go to (71.86) if Parity = 2 


490 


(71.810 93(71.86) JL2F 


From 4 (71 086) 


5 5 




(71.86) KL2F 552042F 


From 10(71.82) 


n 2 : = n 2" 2 




2S1426F 551530F 








2S1530F 854 (j) 








8o4(jl) 8s4(j) 




4(J): = 4(J) -1 




82(71.84) 5F56F 




n^: = n 2 




8KL023F 8I+98OF 








8s4(j) 851004F 




4(j): =^-1 




82(71.88) 93(71.90 


) 


Go to (71° 90) if Parity = 2 




(71.88) 5L2F KL2F 


From 2(71.92) 




500 


(71.90) 511426F 2S1426F 

851004F 83(71.92) 
551530F 92(71.92) 


From 7(71.86) 






(71.92) 511530F 2S1530F 


From l(fl.90) 

» 






8KL023F 844 (J ) 








8s4(j) 83(71.88) 




4(J): = 4(j)-l 




8kf 8^976f 








8s4(j) 2KF 




4(J): = i 




2F57F 5KF 




\: = (0, n 3 ) 




(71.91+) 5F61F 2516 34f 


From 4(71.94) 


ly = (0, n ? ) 


510 


5S1634F 8KL023F 








844 (J ) 8s4(j) 




4(J); = 4(J) -1 


512 


KL2F JL4F 




n 2 : = n 2 -2; ry = ry4 



LOCATION 


ORDER 
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513 


(71.96) 


82(71. 9*0 5K2^F 
851004F 82(71.96) 
52(71.96) 8kf 




b 5 : = (-24,0) 




(71.98) 5Sl63^F 5L2F 

53(71.98) 851004F 




Set odd coefficients = if 

Parity = 2 

F: = 1 ■- Parity 






82(72) 8K2F 










8798OF 8kh(ji) 






520 

■ 


(72) 


.92(72.1) 8K2F 


From 2(71.98) 






(72-1) 


8798OF 8S98OF 




t\ = 2 £ +1 if Parity = 2 




(72.5) 


8KF 8KF 

8n6(N1) 83(72.6) 


From 3(71.82) 


Go to (72.6) if Mode I 




(72,52) 


2K51F 21l63^F 
2Sl63^F 22(72.52) 




e : = -e 
r r 




(72.6) 


8j(20^) 8KL022F 
8^980F 80978F 


From 7(71.70), 
l(71.78),l(72.5 


Go to playback Program A, Ten 
late problem 

BEGIN Elimination of all but 
k+2 Extrema 




(73) 


8S1022F 5KF 
856 (Jl) 8S1012F 


From (9k) 


1022F: = £-k-2 

153 
min _ : = «5 x 10 

r-1 


530 




8KF 8S1014F 
8K2F 8S^(j) 
851022F 80984F 

83(7*0 92(9*0 




max r _ i: = 
MJ): = 2 

F: = t »k-4 




(7*0 


8S1022F 551530F 
83(78) 5H530F 


From 3(83) 


BEGIN Eliminate 2 consecutive 
Extrema or First and Ladt Ext 




(78) 


8S1016F 551532F 
83(79) 5H532F 




min : = |Extremum-. | 




(79) 


8S1018F 8S1020F 




max : = |Extremum | 

r ' 2 J 


539 




8N1016F 83(80) 




F: = max - min 
r r 
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LOCATION 
5^0 




ORDER 

85IOI6F 8S1018F 

851020F 8S1016F 




max : = lExtremum., ( 
r ' 1 

min : = |Extremum_ 1 
r ' 2 




(80) 


85IOI6F 8N1012F 
8S1020F 82(81) 


From 1(79) 


1020F: = min r - miji r .i 




(81) 


92(82.5) 811020F 
83(82) 92(83) 




Go to (82) if min r _ 1 > min r 




(82) 


85IOI8F 8N101^F 








(82.5) 


82(83) 851016F 
8S1012F 851018F 
8S101^F 85^ (J ) 


From (8l) 


Go to (83) if iaax r > max r __ 1 

min , : = min 
r-1 r 

hex ^ = max^ 


550 




8S1010F 2KF 








(83) 


2F59F 5L2F 
8KLF 8H(j) 


From l(8l), 
(82.5) 


V = (0, n^+2) 






8SMJ) 85980F 
8oM«j) 82(7*0 




k(j): = M<J) +1 

Go to (7*0 if t > ^(J) 




(8*0 


$51^30F 82(85) 








(85) 


811530F 8S1016F 
5515 30F 82(86) 




min : = 1 First Extremum| 
r ' 




(86) 


511530F 8S1018F 

8S1020F 8N1016F 




max i = iLast Extremum| 
'r ' 

F: = max r - min r 


560 




82(87) 851016F 
8S1018F 851020F 




max : = iFirst Extremum| 
r ' 




(87) 


8S1016F 851016F 
8N1012F 8S1020F 
83(88) 92(90) 


From 2(86) 


min : = ILast Extremuml 
r ' 

1020F: = min r - min r _x 




(88) 


811020F 82(89) 




Go to (89) if min > min r 


566 


(89) 


93(91) 851018F 







5 LOCATION 




ORDER 
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567 




8N101*fF 8S1020F 








(90) 


85(91) 2KF 
8K2F 8S1010F 


From 2(87) 


Go to (91) if max > max -. 
v ' r — r-1 


570 


(91) 


851020F 85(92) 
252<MF 2S2042F 
25lte8F 2S1426F 
251552F 92(92.5) 


From (89), 
(90), (95) 


Go to (92) if eliminating 
interior extrema 
Eliminate endpoint extrema 




(92) 


251^50F 2S1^26F 
2520^-6F 2S2042F 


From (91) 


Eliminate interior extrema 




(92.5) 


25155^F 2S1550F 
2I2F 8KLF 
8^1010F 8S1010F 
80980F 82(95) 


From 5(91) 


■j f : = J f +l 

c 


580 


(95) 


95(91) 85980F 

8098^F 8S98OF 
8098J4-F niAf 




X : = L -2 




(9M 


92(75) Qhk(ji) 

85(95) 92(99) 


From 5(75) 






(95) 


85155OF 82(96) 




BEGIN Remove One Endpoint 




(96) 


811550F ^R1550F 
32(98) 2KF 
8KL022F 84980F 




Go to (98) if |E(x 1 )| > 
l E ( x f +5 )l 




(97) 


dsk(j) 251^2 8f 


From (98) 


4(J): = f -2 


590 




2S1426F 252044F 
2S20^2F 251552F 
2S1550F 2L2F 




2 X 

x- n : = x^ 
r-1 r 

x .. i = x 

r-1 r 

E(x r _ 1 ): = E(x r ); n^i = n 2 +2 


595 J 




65MJ) 80MJI) 




F: = MJ)-1 

1 
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594 



(98) 



600 



(99) 



(100) 



(101) 



610 



620 



83(97) 8KL023F 
84980F 8S98OF 
NI2F 8KLF 
8S1022F 2KF 
8KLF 8s4(j) 
2S1756F 5KF 
854 (J) 801022F 
8S1020F 8KF 
8N1020F 83(102) 
251426F 501426F 
271736F 2S1736F 
85978F 844 (jl) 
8o4 (J ) 82(103) 
93(104) 8KLF 
844(j) 8s4(j) 
5L2F 93(101) 

8KF 8n4(ni) 
83(104.4) 844(jl) 
83(104,2) 93 (104.4) 

(104.2) 252042F 271736F 
2S1736F 93(104^4) 

(104.4) 85978F 844(ji) 
801022F 82(105) 

(105) 93(106) 8K1F 
841022F 8S1022F 
2I2F 93(100) 

(106) 8KF 8S1022F 



(102) 



(103) 



(104) 



From 1(96) 



From 1(94) 



From 2(105) 



From 2(103) 



From 2(101) 



From (103) 



From l(l04), 
2(104) 



From (105) 



C : = L -l 

BEGIN compute at f . 

V : = lj - V = ° 
4(j): = l 



it 1 : =0511,-: = 
9 5 



1020F: = j-j^ t 

Go to (102) if (j-j^,) = 



tc': = at ' . ..(x. , -x.) 
j-l v jot' ,j' 



Go to (103) if k+l-J > 

Go to (104) if finished with 

this at f . 

4(j)s = J 4(j) +1 



Go to (104.4) if Parity = 
Go to (104.4) if Parity = 2 



F: = k+l-j , 

at 

Go to (106) if finished with 
all at 1 .*s 

n 2 : = n 2 +2 

BEGIN Compute Coefficients of 
polynomial 



LOCATION 



621 



630 



6kO 



6k7 



ORDER 



5KF 8K1F 
8^978F 8sMj) 
8KL023F 8S1018F 

(107) 85IOI8F 83(108) From 4(108) 
8KLF 92(108) 

(108) 8KL023F 8S1018F From (107) 
561736F 84-1022F 
8S1022F 854 (J) 
80k(jl) 8SMJ) 
5L2F 83(107) 
5KF 8KF 
8S1020F 8K1F 

(109) 84978F 8s4(j) From 4(109) 
5515 30F 561736F 

/ 841020F 8S1020F 

854 (J) 8oV(ji) 

5L2F 82(109) 
811020F 861022F 
(109.2) 93(109.5) 50(109.2) From 2(130) 
26 (A7) 8kf 
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(IO9.5) 8S1022F 8KF 
(110) 5K51F 5S1940F 
52(110) 8KLF 
8S1018F 8S1018F 
2KF 8KL023F 
(112) 8S1020F 851020F 

82(113) Self 



From (109.2) 



From (116) 



j : r="k+l 

Binary Switch: = -1 

Binary Switch: = +1 



j: = 3-1 



n^: = n^+2 



a: = 

r: = k+1 or r-1 

a: = a r _ 1 + (E (x r )/it f (x p ) ) 

F: = r-1 

n^: = n^+2 

F: = -a/b 

BEGIN Step I of MODE II 

F: = 

h: = or - a/b 

1940F + ru: =0 
5 

BEGIN store y(x ) + h at 1736 + n^ 
r: = 1 



n 2 : = 



F: = +1 



LOCATION 


ORDER i 
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6kQ 


(113) 


93(11*0 8KL023F 




F: = -1 






(11*0 


8S1020F 871022F 




Binary Switch: = + 1 




650 




2^15 30F 2S1736F 
85978F 801018F 




1736 + n 2 : = y(x r ) + h 

F: = k-r 






(115) 


83(115) 92(116) 
851018F 8^(J1) 




Go to (ll6) if finished storing 

y(x r ) + h 






8S1018F 2L2F 




r: = r+1; n : = n 2 +2 






(116) 


92(112) 8K1F 


From 3(11*0 


BEGIN compute (^+2-^+1 ^ X j+1+ & 






8S101*l-F 1KF 




- x - ,n 






(117) 


1L2F 8KF 
8S1018F 2KF 
2F55F 5KF 


From 10(118) 


n • = 2r 
j: = 
n 2 : = 2f 




660 


(118) 


251426F 501^26f 
8S1016F 551738F 
501736F 861016F 


From 7(H8) 


-- 








5S1736F 5L2F 
2L2F 8K1F 




1736F + 2 .s = ( y . +1 - yj+1 )/:( Xj 


4-1+ 






841018F 8S1018F 




i ' ■" j+i 








85978F 801018F 












801014F 83(118) 




Go to (118) if k- f -j > 








8KLF 8^10l4F 












8sioi4f 85978F 




/?:=/?+! 




670 




80101^F 83(117) 
851736F lS19 i <-0F 
8KF 8S1018F 




Go to (117) if K - ^> 

cr k+l : " z l 
/: =0 






(119) 


2KF 2F55F 


From (123) 


°2 = n l 




67k 




85978F 801018F 









LOCATION 



ORDER 



675 



680 



69O 



700 



(120) 8S101^F 5H^26f 
2719^0F 241938F 
2S1938F 8KLF 
84l01^F 8S101^F 
2L2F 85978F 
801014F 82(120) 
551736F 1^1938F 
1S1938F 8KLF 
&U018F 8S1018F 
9L2F 5L2F 
80978F 82(123) 

(123) 93(119) 5KF 

(12*0 8K1F 8sM<l) 

551634F 50i9if0F 
5S(N) 8KF 
8n6(N1) 83(12^.5) 
55l63^F 54i9^0F 
5S(N) 5S(N) 

(12^.5) 5L2F 85976F 

80^(J) 82(125) 

(125) 93(125.5) 8K1F 
Qkk(j) 92(124) 

(125.5) 856 (Jl) 806 (KL) 
83(126) 8K2F 
87 (Jl) 8S6(N1) 

(126) 8j(206) OOF 
26759N 
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From 5(120) 



From 1(125) 



From 3(12*0 



From (125) 



J:= *-fi 



j: = j+1 
ii g : = 132+2 

Go to (120) if k-j > 



*y =/ /7 + 2 



fix =/?+! 



From 1(125.5) 



n, : = n,-2; n R : = n R +2 



Go to .(123 ) if finished computing 
polynomial coefficients 



j : = 1 or j+1 



<* +1 e^ A ^ +1 

e. : = e. - Ae. 

Go to (12^.5) if Mode II 



^+1 ^ , A -^ +1 
e. : = e" + Ae. 



Fs = k-j; Go to (125) if finished 



Go to (126) if 5 < .5 x 10 155 

8: = .lxlO 155 

Go to Program B, continue Step II 



LOCATION 



ORDER 



00750K 

(127) 50430F 50750F 
26 (Yl) 005953F 
00320F 22(42.62) 

(128) 50^30F 5075 3F 
26 (Yl) 006595F 
00320F 22(71.51) 

(129) 50^30F 50756F 
26 (Yl) 0062 7^F 
00320F 26626F 

(130) 50^30F 50759F 
26 (Yl) OO6595F 
00320F 22(109.2) 
001(207) 02K 
2616F 93(30o05) 
002(204) 02K 
00320F 26(27.1) 
26762N 



NOTES 
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From (71.51) 



From 4(42. 63) 



By 2(30.10), 

16(71.51) 
From (8) 



From 1(71.16) 



Playback Program B 



Go to (42. 6^1) to print out 

extrema 

Playback "Program D 



Go to (71.51) (or 16(71.51) 
if two few extrema) 
Playback Program C 



Go to (71°3) to read in Auxiliary 

Subroutines 

Playback Program D 



Go to (109.2) for first iteration 
of Mode II . 



Go to (30.05) to continue Step II 
iteration 



Go to (27.1) to terminate problem 



